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SUMMARY 


A theory is initiated for the generation of waves upon a water 
surface, originally at rest, by a random distribution of normal 
pressure associated with the onset of a turbulent wind. Corre- 
lations between air and water motions are neglected and the 
water is assumed to be inviscid, so that the motion of the water, 
starting from rest, is irrotational. It is found that waves develop 
most rapidly by means of a resonance mechanism which occurs 
when a component of the surface pressure d'stribution moves 
at the same speed as the free surface wave with the same wave- 
number. 

The development of the waves is conveniently considered in 
two stages, in which the time elapsed is respectively less or greater 
than the time of development of the pressure fluctuations. An 
expression is given for the wave spectrum in the initial stage of 
development (§3.2), and it is shown that the most prominent 
waves are ripples of wavelength A,, = 1-7 cm, corresponding to 

. the minimum phase velocity c = (4g7/p)'! and moving in 
directions cos-(c/U,) to that of the mean wind, where U, is 
the ‘convection velocity’ of the surface pressure fluctuations of 
length scale 4, or approximately the mean wind speed at a height 
., above the surface. Observations by Roll (1951) have shown 
the existence under appropriate conditions, of waves qualitatively 
similar to those predicted by the theory. 
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Most of the growth of gravity waves occurs in the second, or 
principal stage of development, which continues until the waves 
grow so high that non-linear effects become important. An 
expression for the wave spectrum is derived (§ 4.1), from which 
follows the result 


2v2p?U. g’ 
where €* is the mean square surface displacement, P? the mean 
square turbulent pressure on the water surface, ¢ the elapsed 
time, U,, the convection speed of the surface pressure fluctuations, 
and p the water density. ‘This prediction is consistent with 
published oceanographic measurements (§ 4.3). 

It is suggested that this resonance mechanism is more effective 
than those suggested by Jeffreys (1924, 1925) and Eckart (1953), 
and may provide the principal means whereby energy is transferred 
from the wind to the waves. 


1. INTRODUCTION 

In the following pages, a theory is proposed to account for the generation 
of waves by wind. ‘The problem has attracted attention for many years, 
during which time many experimental studies have been made and many 
hypotheses advanced to account for the results of these observations. 
However, no really satisfactory explanation of the phenomenon has been 
offered and not even the physical processes involved can be regarded as 
known. Although casual observations of wind-generated waves must have 
been made by almost everyone, thorough experimental investigations are 
difficult, and our understanding of the subject has consequently been very 
limited. ‘The best known of the various hypotheses put forward are those 
due to Jeffreys (1924, 1925) and Sverdrup & Munk (1947), which are 
described in a valuable review of the subject by Ursell (1956).  Jeffreys’s 
sheltering hypothesis is the very reasonable one that the wind passing over 
any waves already existing on the surface induces a variable pressure 
distribution, and that the component which is in phase with the wave 
slope supplies further energy to the wave, resulting in its continued develop- 
ment. ‘lhe predictions from this hypothesis involve a sheltering coefficient s, 
which, for a wave in which the energy gain from the wind just balances the 
viscous dissipation, is required to have a value of approximately 0-27. 
However, experimental measurements of s made on solid wave models 
gave values which were smaller by an order of magnitude, so that it must 
be concluded that the etfect of sheltering, though undoubtedly present to 
some degree, is not efficient enough to account for the observed wave growth. 
The failure of this hypothesis led Sverdrup & Munk and others to consider 
the etfect of the tangential stresses set up on the water surface. By assuming 
that a// the energy communicated to the water by these stresses appears in 
waves and none in currents, they were able to account for the order of 
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magnitude of the observed wave heights. But such theories, based upon 
energy considerations alone, are hardly satisfactory; and, as a means of 
producing the irrotational motion associated with surface waves, the shear 
stresses seem intuitively to be singularly inefficient, particularly in the 
early stages of wave growth. It seems more plausible to suppose that the 
initiation and early development of waves is a consequence of fluctuations 
in the normal pressure upon the surface, which under natural conditions 
are always present on account of the invariably turbulent character of the 
wind. As the waves develop, the influence of the shear stresses upon the 
wave pattern may not remain negligible, though it is difficult to give an 
a priori estimate of their importance. However, it will be seen later that 
the observed wave heights can be explained in terms of the pressure fluctua- 
tions alone, so that it seems possible that the effects cf the shear stresses 
are rarely, under natural conditions, dominant in the wave growth. 

The present approach is developed without specific empirical assumptions, 
such as the sheltering hypothesis of Jeffreys or the assumptions of energy 
transfer involved in the work of Sverdrup & Munk. ‘The mathematical 
problem can be formulated as follows. Given that at an initial instant a 
turbulent wind commences to blow across an infinite sheet of deep (inviscid) 
water originally at rest, generating a distribution of fluctuating pressure on 
the surface which is a stationary random function of position and time, the 
aim is to study the properties of the surface displacement at subsequent 
times. 

An attempt to solve the problem in a manner similar to this was made 
by Eckart (1953). He represented the pressure distribution associated with 
the turbulent wind in a very specific way, i.e. as an aggregate of pressure 
points of a given size and duration over a finite storm area. It is very likely 
that this representation is too specific to give an adequate description of the 
randomness of the fluctuations in pressure associated with the turbulent 
wind. It will be seen later that this randomness is essential, since the water 
surface selects from the whole range of wave-numbers in the applied pressure 
distribution certain ones whose amplitude increases most rapidly by a 
mechanism involving a type of resonance. ‘The inadequacy of Eckart’s 
representation in this respect is probably the most important cause of the 
failure of his theory to explain the observed wave heights; his predicted 
values are an order of magnitude too small. 

It is clear from the way the problem has been formulated that the present 
treatment is not a stability theory, such as the Kelvin-Helmholtz model 
(see, for example, Ursell 1956). It is found that waves can grow not only 
through an instability but by this mechanism of resonance between the 
components of the surface pressure distribution and the free surface waves. 
It cannot be claimed that this resonance effect, whose nature will be 
explained more fully below, is the only one capable of generating waves; 
for instabilities of the Kelvin-Helmholtz type will almost certainly occur 
if the wind velocity is sufficiently high and established sufficiently rapidly 
(perhaps producing the capillary ripples observed in strong gusts of wind), 
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and in some circumstances the effects of sheltering and shear-stress action 
may not be negligible. However, it does account in a natural way for a 
variety of observed phenomena, though a true appreciation of its scope and 
limitations must await the results of further investigations, particularly on 
the effects of the non-linearities which become important when the wave 
slopes become appreciable. 

My understanding of the problems involved in the generation of waves 
by wind has been greatly helped by discussions with Dr M. S. Longuet- 
Higgins of the National Institute of Oceanography, Surrey, and with 
Mr F. Ursell and Dr A. A. Townsend at Cambridge; it is a pleasure to 
acknowledge my indebtedness to them for their constructive comments 
upon this work. 


2. ‘THE APPROACH TO THE PROBLEM 
2.1. The physical background 

It is a matter of common observation that winds, blowing over the 
ground or over water, do not consist of streams of air in steady and uniform 
motion but rather of an irregular series of ‘ puffs’ and ‘lulls’ carrying eddies 
and swirls distributed in a disordered manner. In the language of the fluid 
dynamicist, the wind is not laminar, but invariably turbulent. This 
turbulent nature of the wind will be found to result in the birth and growth 
of waves upon a water surface. The atmospheric eddies, or random velocity 
fluctuations in the air, are associated with random stress fluctuations on 
the surface, both pressures (i.e. normal stresses) and tangential shear 
stresses. ‘The eddies are borne forward by the mean velocity of the wind, 
and at the same time they develop, interact and decay, so that the associated 
stress distribution f both moves across the surface with a certain convection 
velocity dependent upon the velocity of the wind and evolves in time as 
it does so. An overall convection velocity of the stress fluctuations can be 
defined mathematically as the velocity of the frame of reference in which their 
frequency scale is least, or their time scale greatest. This least frequency 
scale is characteristic of the evolution of the stress distribution, whereas 
the frequencies observed, say at a fixed point, are invariably greater because 
of added contributions from the sweeping of the stress pattern past the 
point. 

The stress distribution contains components with a very large range of 
wave-numbers, and as might be expected, Fourier components with very 
different wave-numbers are convected at different speeds. ‘The surface 
stress fluctuations of any given length scale can be attributed to turbulent 
eddies near the surface whose length scale is of the same order of magnitude, 
and the larger eddies will be convected by the mean wind with a greater 
speed than the smaller stress-generating eddies lying near the surface. 
A definition to express this slightly more refined concept can be given as 
follows. If F(x, ?t) is the two-dimensional spatial Fourier transform of the 


surface stress covariance f(x’, /’)f(x’ + x, +1), then the convection velocity 
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U(x) for the components of wave-number x is the velocity of the frame 
of reference in which the time scale 


[ F(x, 0)? | F(x, t) dt 
“0 
is greatest. ‘l‘hat is, we take axes moving with velocity u, and in this frame 
of reference (see § 4.1) 


F(x, t,u)=[F (x, t)]y — cos(x. ul). 


The convection velocity U,(%) is the value of u tor which 

| — cos(x.ut) df 

is a maximum. ‘l‘his maximum time scale, or development time of the 
stress fluctuations, which appears later in the analysis, will be denoted 
by ©(x). In the frame of reference convected with velocity U,, the time 
dependence is slowest and is determined by the growth and decay of the 
eddies as they are carried along. 

The convection velocity U, defined in this way must be related to the 
mean velocity of the wind measured at a certain height above the surface. 
‘The physical picture of the convection of the stress-generating eddies 
mentioned above suggests that U.(x) is of the same order as the mean 
velocity of the wind at a height «' above the surface. Now, the mean 
velocity profile of the wind U(z) is determined by the friction velocity 
ux (=(7»/p,)'” where 7, is the mean horizontal shear stress at the surface, 
and p, the density of the air), by the roughness length z,, and by the 
kinematic viscosity v, and is of the form 


U(z) 


Uy 


In aerodynamically rough flow, the function B is in effect zero, and the 
roughness length 2, for a water surface is probably related to us by an 
expression of the type uZ & gz) and for moderate winds is of order 0-1.cm. 
The velocity profile of the wind over a water surface is discussed in recent 
review articles by Ellison (1956) and Ursell (1956); the reader is referred 
to these accounts if further detail is required. However, the important 
point in this connection is that when 2/2, is large, the variation of mean 
velocity with height is slow, so that for the components of the stress 
distribution whose length scale is of the order of several metres, the 
convection velocity U(x) is very nearly equal to the wind velocity as usually 
measured from a ship. For the very much smaller-components, however, 
it may be considerably less. 

Our problem is to discover the nature of the action between this stress 
distribution and the water surface, and the investigation of this is described 
in the following sections. However, it will perhaps be clearer if some of the 
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salient features are described first in physical terms, although this will 
necessitate anticipation of some of the analysis to be developed later. ‘This 
background will, it is hoped, make the direction and implications of the 
subsequent mathematics more evident. 

We have assumed that the fluctuating pressure upon the water surface 
is responsible for the birth and early growth of waves. If a Fourier analysis 
were to be made in the two space dimensions of the surface and in time, 
components would be found over a wide range of wave-numbers and 
frequencies. ‘hese components of the pressure fluctuations acting upon 
the surtace generate small forced oscillations, and the amplitude of any 
Fourier component of the surface displacement depends upon the amplitude 
of the corresponding component of the forcing pressure fluctuations. 
But the response of the water surface to the various Fourier components 
of the pressure field is not uniform, since certain combinations of wave- 
number and frequency in the components of the surface displacement are 
excited more readily than others. In particular, if the pressure distribution 
contains Components whose wave-numbers and frequencies coincide with 
possible modes of free surface waves, then a type of resonance occurs, and 
the continued presence of these components in the pressure distribution 
(although their amplitudes may not be large) generates surface waves whose 
amplitude continually increases. 

The nature of this resonance can be seen most clearly by considering 
the motion in a frame of reference moving with the convection velocity 
of the surface pressure distribution. Suppose we can, for the moment, 
neglect the evolution of the pattern of stress fluctuations so that the 
distribution is stationary in the convected frame of reference. ‘The various 
components generate wavelets of small amplitude which propagate along 
the surface in all directions; but for a given convection velocity (provided 
it is not too small) there are two particular wave modes that propagate in 
the direction of the wind at the convection velocity. These particular 
wavelets remain stationary with respect to the pressure distribution which 
produced them. ‘lhe phase ditference between these particular wavelets 
and the Fourier component of the pressure field with the same wave-number 
remains fixed, so that the conditions which gave them birth remain to 
ensure their continued growth. 

But this is an over-simplification, because the evolution of the stress 
pattern usually cannot be ignored. If we make a two-dimensional spatial 
Fourier analysis of the surtace pressure distribution with respect to axes 
moving with the convection velocity, the amplitude and phase of each 
component varies slowly in a random manner with time. ‘The particular 
wavelet which is stationary with respect to the convected reference frame 
is still subject to a pressure component of the same wave-number, but one 
whose phase now wanders and whose amplitude varies. However, the 
average amplitude of the wavelet continues to increase, although we should 
expect the variation of mean square amplitude at this wave-number to be 
slower than it would be if the phase were constant. ‘This process of 
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statistical development through a random influence is already familiar in 
the contexts of diffusion and random-walk problems. 

The mathematical problem in which this process appears has already 
been formulated in the previous section, and certain assumptions will now 
be made to facilitate the analysis. We should like to be able to assume 
that the statistical properties of the applied pressure fluctuations are 
independent of the waves already generated. We will therefore endeavour 
to determine the conditions under which this is likely to be true. Initially 
the water is at rest, so that in the very early stages of wave growth (i.e. the 
‘initial stage’ of § 3) this assumption can clearly be made as a first approxi- 
mation, since the initial pattern of ripples has not yet grown large enough 
in amplitude to have a significant effect on the wind. Presumably, however, 
the pattern of ripples soon reaches some sort of statistical equilibrium, and 
there is evidence (see Ursell 1956, p. 240) that the roughness length 2, 
of the surface is determined by these short steep ripples rather than by 
the height of the longer waves that subsequently develop. It seems, then, 
that after quite a short duration these ripples begin to exercise a profound 
influence upon the mean velocity profile and the other mean properties 
of the turbulent wind. However, a much greater time is taken for the 
development of the longer gravity waves in which we are more interested. 
It seems, therefore, that when the duration is sufficiently long, the 
components of the pressure distribution generating the longer gravity 
waves are independent of the waves of this kind already present. ‘To 
summarize, this assumption seems justified under two sets of conditions: 
(i) when the duration of the wind is very small, and (ii) for gravity waves 
when the duration is large. In both cases, this assumption must be 
expected to fail when the wave steepness becomes too great, since, for 
example, sheltering effects may then become important. 

As a further simplification, the viscosity of the water will be neglected, 
since it 's probably unimportant for all but the shortest waves over moderate 
intervals of time. The motion, which is envisaged as starting from rest, 
is then irrotational. Consideration is restricted to waves whose mean- 
square slope is small, so that the surface boundary condition can be linearized. 
Finally, for definiteness and convenience, the surface pressure fluctuations 
are assumed to be statistically homogeneous over regions of the surface 
whose extent is comparable with the correlation areas of the atmospheric 
turbulence, and, after the initial instant, to be statistically steady. 


2.2. The equations of motion 

In this problem it is convenient to express the surface pressure fluctuations 
and surface displacements in terms of their Fourier components. The 
pressure fluctuation p and displacement € are stationary random functions 
of position x = (x, x.) in the surface plane, so that simple Fourier transforms 
do not exist. However, we can invoke a more effective technique, used 
widely in the theory of random functions (see, for example, its application 
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to turbulence theory in Batchelor (1953)), by defining the Fourier—Stieltjes 
transform 


E(x,t) = | Nem, (2.2) 


where the integration is over all wave-numbers x in the plane. The two- 
dimensional instantaneous spectrum of the surface displacement is the 
Fourier transform of the covariance 


= €(x)&(x+r), 


which decreases to zero rapidly as r -> ~, so that its Fourier transform 


(x) = (27)? | dr, 


where dr represents dr, dr, and the integration is over the entire surface 
with x fixed, does exist in the ordinary sense. From the inverse of this 
relation, it can be shown that the spectrum is given in terms of the 
Fourier—Stieltjes transform by 


W(x, t) = dA(x, t) dA*(x, “J (2.3) 


dk, dks 


where the asterisk indicates the complex conjugate. Similar expressions 
can be written down to specify the pressure fluctuations on the surface. 
The two-dimensional Fourier—Stieltjes transform 


p(x, t) = | (2.4) 


is related to the spectrum function II(x,¢) of the surface pressure 
fluctuations by the expression 


2) = (22)-* | p(x, t’)p(x+r, t' +t) e™-* dr, 


dex(x,t) da*(x,t +1) 


dk, 


wm 


which is a function of the wave-number x and time separation ¢. 
The motion in the water has been assumed to be irrotational, so that 
a potential ¢ can be defined such that 


and V2¢=0, (2.6) 


where u is the velocity vector. ‘The surface boundary condition for waves 
of infinitesimal height is 


OF od 


where p is the density of the water, 7 the surface tension at the interface, 
g the gravitational acceleration, € the surface displacement, and z is the 
vertical position coordinate to be measured downwards from the undisturbed 
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surface level. It is convenient to express this equation in a frame of reference 
moving with an arbitrary velocity U, and it becomes 
- =] —-U,— —gé+ + (2.8) 
p ct = 6 p 
The normal velocity of the surface is d&/dt, or in the moving frame of 
reference, 


Ux) B | (dA’—ix.U dA)e**, (2.9) 


Ca 
where the representation (2.1) is used, and where the accent represents 
time differentiation in the convected frame of reference. ‘The irrotational 
motion in the water is specified by the solution of V*¢ = 0 with the boundary 
conditions that dé - 0 as s—> %, and that at s = 0, 0d/dz is as given by 
(2.9). ‘This solution is (see Phillips 1955 

dd’ —in. Uda | 


K 


(2.10) 


where « = x. ‘The equation (2.8) can now be expressed in terms of the 
Fourier—Stieltjes transforms as 


p ct kK p 
or dA” —2in, dA’ —(n? —n3) dA = — de), (2.12) 
where n, =xU cosa, Ry = (ge + Tk?/p)'?. (2.13) 


Equation (2.12) describes the growth of each component of the surface 
displacement in terms of the corresponding component of the pressure 
distribution. If we now specify U as the convection velocity of the component 
of the pressure field of wave-number x, the quantity m, represents the 
frequency, in radians per second, of a wave with wave-number « and speed 
U.x/x = U,.cosx in a direction at an angle « to that of the wind, and a, is 
the frequency of free surface waves of wavelength 27/k. 

The fundamental equation (2.12) specifies the development with time 
of the Fourier—Stieltjes transform dA(x, t) of the surface displacement in 
terms of the transform da(x,?t) of the pressure fluctuations, which is, 
after the initial instant, a stationary random function of time. Its solution, 
subject to the initial conditions appropriate to an undisturbed surface at 
rest, namely dA = dA’ = 0 at t = 0, can be expressed as 


ells rt K 
or 

dA deo(s)[exp{ —i 
A(x, t) = 2pn, w(7)[exp{ —n2)(7 —t) — 


dr. (2.14) 
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This solution provides the common starting point for the investigations 
of the next two parts of this paper. It illustrates how the amplitude of each 
Fourier-Stieltjes component of the wave pattern is dependent upon the 
history of the applied pressure distribution from the initial instant t = 0. 
The development of each component of the wave field falls naturally into 
two stages: the intial stage over values of time t < (x) (i.e. the development 
time of the fluctuations in pressure of wave-number x, or the time scale 
of dw(x, t) in the convected frame of reference), and the second or principal 
stage of development for which t > © and during which most of the growth 
of gravity waves takes place. ‘The properties of the waves in each of these 
stages will be discussed in turn. 


THE INITIAL STAGE OF WAVE DEVELOPMENT 
3.1. The initial response of the water surface 
We define the initial stage of development as the period of time for 
which t < @(x), so that over this range of t, do(x,t) in the convected 
frame of reference varies little from its value at the initial instant. In the 
integration of the general solution (2.14), dw(x) can be considered indepen- 
dent of time, and the solution becomes 


ik du(x) 
p(n? — n5) | exp{ —1(m, + 
— Ny 


From this, there follows a relation between the wave spectral function @(x) 
and the two-dimensional, simultaneous pressure spectrum I1(%) = I1(x, 0); 
for 
dar*(x) nN, +My Ny — Ny 
= 2 exp!i(n, — 
<< | exp{i(m, — + =eXxp{i(m, +ny)t} 


— 2n, 


dA*(x, t) 


and since 


_ dA(x)dA*(x) du(x) da*(x) 
we have 
(D(x, t) == — + = —Ny)t + = cos(m, + — 


~ F 
cos 2n, (3.2) 


If we write 


vy, = n,t = Vy = Nyt = (ge + Tk?/p)!? (3.3) 
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a dimensionless response factor can be defined as 


= (v7 — v3) *4 cos(v, — cos(v, + Vy) — 
2 2 
ve 
(3.4) 


so that the wave spectrum can be expressed in the form 

D(x, t) == (x)t/p?. (3.5) 
The tunction (x, 7) is therefore determined, when ¢ < ©(x), jointly by 
the spectrum of the imposed pressure fluctuations I(x) and by the response 
factor I’, itself a function of « through v, and v,, and of t; the balance 
between these two factors must always be taken into account in considering 
the wave generation in the initial stage. The function «?II(x) has its 
maximum when « is large, of the same order as the wave-numbers associated 
with the dissipating eddies in the atmospheric turbulence. A brief examina- 
tion of (3.4) suffices to show that in the (v,, v))-plane, [' is large near the 
line v,; = v, and in regions where v/v, 1s large. 

Wave-numbers of the first class, for which v, = vs, are such that the 
convection velocity of the pressure distribution is equal to the velocity 
of free waves with these wave-numbers in the direction of convection ; 
tor, from (3.3), 

U..cos% = (g/k+ Tk/p)!? = e(«). (3.6) 
Because of the factor «*II(x), the amplitude of these waves is significant 
only when « is large (i.e. the wavelength A is small), which in turn implies, 
from (3.6), that U’.cosx should also be small. Waves, or perhaps more 
properly ripples, of this type are generated immediately after the onset 
of a gust of sufficient strength. ‘lheir direction and speed is such that they 
remain in phase with the applied pressure distribution during this initial 
period, developing by the simple resonance mechanism described in § 2.1. 
‘These waves can conveniently be called ‘ resonance waves’. Wave-numbers 
of the second class, for which v, > vy, are associated with the forced 
oscillations of the water surface produced by the pressure fluctuations, 
and it will be shown later that their amplitude does not increase with 
time ¢ beyond a small upper bound. 


Resonance waves 
A simple approximation to the response factor (3.4) when », = vy 
can be found in the following manner. Let 
X= ¥y— (3.7) 
and, on substituting this into (3.4), we have 
45 vet 2vex+ x” 


= (2v,+ x) cos x + * cos (2v,+x)- 
2 dv, A 
_ 


2 


as 
— COS ¥ \ } 
Vo 
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If v, > 1, that is, if t > n>! or (27)! times the period of free waves of 
wave-number «, the second term is negligible and the wave spectrum is 
given approximately by 

1—cos x 


D(x, t) = (3.8) 


when ny! <t < © and n, = ny. 

‘The existence of a range of values of ¢ satisfying this condition clearly 
requires that ny! is less by one or two orders of magnitude than 0. Since 
n, — My, tor resonance waves, an equivalent statement is that the wave- 
length A of the surface disturbance is very much less than the distance 
travelled by the convected pressure distribution in the time ©. ‘This 
condition is likely to be satisfied best when A is small, in the range corre- 
sponding to ripples and short gravity waves. (If, for example, (x) is 
of the same order as the time scale of the eddies of wave-number « in the 
equilibrium range, then the condition can be expressed as 


A= 2alk < U.O Ue 


(where ¢ is the rate of energy dissipation) or « > 87%e/ U7.) 
‘The wave-number « in any direction « for which ®(x) is a maximum 
is determined by the equation y = 0, or, from (3.3), 


Tx?/p —« Uz cos?x +g = 0, (3.9) 


which, for each value of U.cosx above a certain minimum, gives two 
resonant wave-numbers, one corresponding to a gravity wave and one 
to a capillary wave. ‘This is illustrated in figure 1. The solid curve c(«) 
gives the speed of free surface waves as a function of wave-number « and 
the broken line U,(«) a typical convection speed of the components of the 
pressure field of wave-number «, which, as we have seen earlier, is a function 
which may decrease slowly as x increases. ‘The exact shape of U,(«) depends 
upon the surface roughness, and the curve in figure 1 is for illustrative 
purposes only. In the direction x = 0, that of the convection velocity, 
resonance waves occur near the wave-numbers such that U,(«) = c(«), 
given in figure 1 by the points 4 and B. As « increases, since the condition 
for resonance waves is that U.(«)cos x = ¢c(x), the resonance wave-numbers 
are found by translating the dashed curve downwards, and the intersections 
of the two curves approach each other until at a certain critical angle x,, they 
coincide. For values of « greater than ~,,, no resonance waves are possible. 
If, as one might expect under normal atmospheric conditions and with az 
‘aerodynamically-smooth surface’ appropriate to the initial period of 
development, the dependence of U, upon « is not very rapid except when « 
is large (of order unity), the coincidence point will occur near or a little 
to the right of G in this diagram, which marks the minimum velocity of 
free waves under the influence of both gravity and surface tension. It can 
readily be established from (3.6) that this minimum velocity is given by 


Cmin = p)! (3.10) 


- 
i 
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so that the wavelength of the critical waves is approximately 
Nee = 20(T/pg)"?. (3.11) 
The critical angle ~,,, i.e. the maximum angle from the direction of the wind 
at which resonance waves are generated, is given by 
Ler = COS "'(Cmin/U,) = (3.12) 
where U, is the convection velocity of the wave-numbers «, i.e. the mean 
wind velocity at a height of approximately «'. ‘Taking the values for 
water: p= 1gmcem~*, 7 = 73 gmsec”, and taking g = 980cmsec?, we 


have Cmin = 23cmsec™! and the critical wavelength 
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Figure 1. 


Of these possible resonance waves, not all may be excited, since for 
some wave-numbers «7II(%) may be negligible. It is well known that in 
turbulent flow there is a certain maximum wave-number kmax of the velocity 
fluctuations, or minimum eddy size, which depends upon the rate of kinetic 
energy dissipation per unit volume in the flow near the surface and upon 
the kinematic viscosity v of the fluid. If the mean velocity gradient is not 
too large, Kmax ~ (e€ v*)!4. The surface pressure fluctuations will also 
have a minimum scale, or maximum wave-number, of the same order of 
magnitude, so that for values of « greater than kmax, II(x) is negligible. 
Not a great deal of relevant experimental information is available con- 
cerning the magnitude of « under the conditions with which we are 
concerned, but for a mean wind speed of order 1 msec! at a height of 
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1 metre, a rough estimate guided by experiments such as those of ‘l'aylor 
(1952) suggests that the maximum value of « is of the order of 10cm"! 
or less. ‘Thus, in figure 1, only those wave-numbers below a band such 
as DE are excited. In the direction x = 0, running before the wind, the 
only resonance waves that occur are the gravity waves corresponding to 
the point B; the wavelength of the capillary waves given by the point 4 
is too small tor excitation to occur. Even if the turbulence were more 
vigorous, with a smaller minimum eddy size, the growth of the very short 
capillary waves is hindered by direct viscous dissipation which becomes 
important for very short wavelengths. Under the conditions illustrated 
in this figure, no true capillary waves at all would be generated, the smallest 
wavelengths being little less than 4,,,. 

It will now be shown that the angle z,, is associated with limiting 
behaviour in yet another respect. Not only does it represent the maximum 
angle from the wind direction at which resonance waves can be generated, 
but also it is found that the amplitude of waves moving in this direction 
is greater than those in other directions with « < z,,. ‘The curve x = 0 
gives the locus of the points in the wave-number plane for which the 
maximum amplification of resonance waves can occur, and for a given 
value of LU’, is of the form shown in figure 2, where the 1-direction is taken 
as that of the wind. For angles z > z,, there are no resonance waves, 
and for angles x — %,, there are two possible values of « for which resonance 
may occur, although one may not be excited if the minimum scale of the 
turbulence is too large. ‘(he maximum wave-numbers that occur in the 
surface pressure distribution can be represented by the are AB, and only 
the wave-numbers near the curve y(x) = 0 and within the sector OAB 
(and its mirror image in OA) will be excited as resonance waves. The 
width of the resonance band is, as (3.8) shows, determined by the range 
of values of y for which (1 —cos x)? is significant. ‘The limits of this 
range can conveniently be taken as y(x) = +27, where, from (3.7), 

= tik U.cos %— (gk + Tx? /p)'?}, (3.13) 
and these limits, when #,¢ > 1, are illustrated by the dotted lines in 
figure 2. It can be shown from (3.13) that the thickness A« of the resonance 
band at the wave-number « is such that Ax <« when ¢t > ny! and 
decreases as t!. As time goes on, the range of excited wave-numbers 
near the curve y(%) = 0 becomes narrower and narrower while the magnitude 
of the wave spectral function (x) along this curve rapidly increases. 

Now the mean-square surface displacement is given by 


= | O(n) de, (3.14) 


the integration being over the whole x-plane. But the integrand is 
significant only in the part of the strip near the curve x(x) = 0 lying within 
the quadrant OAB. A simple and useful approximation to the integral 
can be obtained by taking locally orthogonal coordinates s, being the 
distance along the curve x(x) = 0 from the point D, and y specifying the 
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distance normal to the curve. 


Since the effective width of the resonance 
band is small, the Jacobian in the expression 


5) 


dyds (3.15) 


Ke 
Figure 2. 


can be evaluated along y = 0. The limits of integration of x can be taken 
as + «x. From (3.13), it can be shown that, when x = 0, 
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The contribution to the mean-square surface displacement per unit length 
of the resonance band is simply 


K, t ie 1—cos 
4 0 


(x, 5) p?U.Dcose J _ « 


27 
and the contribution per unit angle in the wave-number plane, or the 
directional distribution of wave amplitude in the physical plane, is found 
by multiplying this expression by es/¢x. With the aid of (3.16), this gives 
2m (x).t 


p?U.cos — Tk?/p)’ 


W(x) = 


where « and « are related by «U,,cos x = (gx + Tx/p)!”. 

From this result, (3.17), most of the important properties of resonance 
waves in the initial stage can be found immediately. In directions z < ~,, 
the resonant wave-number «x < x,, = (pg/T)!*, so that (3.17) becomes 
approximately 
I1(x)t 


% 


= (3.18) 
However, as « —--«,,, T?/p > g, and from (3.17), the directional distribution 
function ‘I'(x) becomes infinite at x... When («,,—«) is small, (3.17) gives 
mK 11(x)t p 
But from (3.13), on the curve x = 0 near ~,,, 


r\P Ce, 


sin 2%, g 


4, the directional distribution function 
of the wave motion is given approximately by 


\ gp sin?2x,, 


so that near the critical direction z 


The form of ‘f(x) is shown in figure 3. As x increases from zero towards 
x... ‘I'(x) increases towards its (integrable) infinity; and when x > ~,,, 
(a) 1s zero to this approximation. 

‘The mean-square surface displacement is obtained by integrating ‘’(«) 
over all possible propagation angles x and wave-numbers « subject to the 
condition y = 0. Nearx = x,, the directionai density of the wave distribution 
is large and should produce wave patterns more pronounced in this direction 
than in others. The reason for this accumulation of waves in the direction 
x,, is, of course, that the range of wave-numbers that propagate at speeds 
within, say, dc of the minimum speed of surface waves is proportional to 
(5c)'? only; within an angle dz of ~,, the range of wave-numbers propagated 
is proportional to (dx)'?, so that the directional density of the wave pattern 
is proportional to 
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The angle between the direction of propagation of these waves and 
that of the wind, i.e. %er = cos~(¢min/U,), will be small only for light airs. 
It is doubtful whether a wind of, say, 30 or 40cmsec™! would contain 
turbulent components whose scale is as small as the 1-7 cm necessary to 
excite such waves, and unlikely that they would be observed at small angles 
%er under natural conditions. As the wind velocity becomes larger and 
contains smaller-scale turbulent components, xr also increases, and these 
critical angle waves may travel almost at right angles to the wind! Under 
the conditions visualized in this paper, with a large fetch and sudden onset 
of the wind, this wave pattern would be a transient one, and would 
disappear as the waves developed beyond the initial period. 


PROPORTIONAL 
TO (Xo - x ) “3 


Figure 3. 


However, it is likely that the waves are similar to those generated in 
steady winds at a small fetch, since the time of development of each wave 
is small, and the growth would be analogous to that discussed here. 
Assuming a correspondence between the two cases of short wind duration 
and short fetch, an explanation of some of Roll’s (1951) observations can 
be given. These experiments were made under natural conditions on 
pools left in tidal mud flats at Neuwerk. On one occasion, when the wind 
velocity at a height of 35 cm was 195cmsec~!, Roll observed a rhombic 
wave pattern of wavelength about 1-7 cm. Unfoftunately he does not 
record the shape of the rhombic pattern, but this observation seems 
consistent with the establishment of the two trains of resonance waves at 
angles x, to the wind direction, as predicted by the theory. On two other 
occasions, when the wind velocity was approximately 40cm sec~', waves 
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with A =1-:7. cm were observed with their crests nearly parallel to the 
direction of the wind, which again is consistent with the theoretical 
predictions. However, many more careful experiments will have to be 
made before quantitative comparison with this part of the theory will be 
possible. 


3.3. The question of Umin 

It is perhaps appropriate at this point to offer some remarks concerning 
the idea that there exists a minimum wind speed capable of raising waves. 
It is natural that this concept should arise if the problem is considered as 
one of stability and many attempts to determine Umin experimentally 
have been made in the past. Ursell (1956), however, discounts the idea, 
pointing out that any turbulent air motion will necessarily generate some 
surface displacements, however minute, and that there seems to be no 
reason why this should not be described as a wave motion. He quotes a 
table of ‘measured values of Umin’ ranging from 40 cmsec~! (Roll 1951) 
to above 1200 cm sec"! (Keulegan 1951). 

According to the present theory, the minimum wind velocity capable 
of raising resonance waves is 23 cmsec ', but waves will actually be generated 
only if there exist near the surface turbulent fluctuations of sufficiently 
small scale. But, as Ursell points out, turbulent motion in the air will 
generate surface displacements irrespective of the mean wind velocity. 


For example, if U). = 0, then in (3.4) v, = 0 and 
1 —cos v,\2 
(3.20) 
v5 
so that, when f < ©), 
(P(x, t) = — cos 
1 
and | ome —cosn,t)? dx 
ny 
4 (x 
5 ii dx = const. (3.21) 


‘The mean-square displacement of these forced waves is therefore bounded 
and unlike the resonance waves they do not continue to develop in time. 
To the eye they may appear as a shimmer or ruffle on an otherwise glassy 
surface. When v, > vy, a similar system of forced oscillations can also be 
found. 


4. "THE PRINCIPAL STAGE OF DEVELOPMENT 
4.1. The wave spectral function 
The principal stage of development for a wave-number x is defined as 
the period during which the time elapsed from the initial instant is greater 
than (x), the development time of da(x,t), but not so great that the 
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amplitude of the wave at this wave-number has increased to such an extent 
that the infinitesimal wave theory ceases to apply. 
We return to the fundamental solution (2.14), namely 


dA(x, t) = der(x, 7)[exp{ —i(n, —n,)(7 — 
2/70 
—exp{—i(m,+n.)(7—t)}] dr. (4.1) 
The wave spectral function ®(x, ¢) is given by 


dA(x, t) dA*(x, t) 
dk, dk, 


D(x, t) 


+exp{—i(n, —7’)} — 2 exp{ —in,(7 —7’) }exp{ —in,(7 +7’)} x 
x cos 2n,t}drdr’. (4.2) 


If the variables of integration are changed from 7,7’ to 7,, 7), where 
7, =7-7', T2=T+7’, it is found that the asymptotic form, as t > 0, 
of the expression (4.2) is 

Kt 


06) ~ 


| II(x, 7){exp{ —1(m, + + exp{ —1(m, — dr, 
(4.3) 


where the suffix has now been dropped from 7,. This asymptotic form 
can be written in a simpler manner if we interpret the integral in terms of 
time scales in frames of reference themselves moving with respect to the 
convected frame in which the analysis has been made. If a spatial Fourier 
transform F(x, 7) of some covariance f(r, 7), a function of space separation r 
and time interval 7, is defined as 


flr,7) = | F(x, de, 


then, in a reference frame moving with a relative velocity V, such that 
r= q+ Vr, we have 


T) = | F(x, r)e™ V2) dx, 


so that 
F(x, 7) = (20)? | dq, 
the q-coordinates being in the moving frame of reference. Thus the 
spectral function in the moving frame is 

F(x, 7), 
where F(x,7) is the spectral function in the original frame. The integral 
time scale, for the wave-number x, in the moving frame, is therefore given 
by 4(x, V), where 


eo 


F(x, 7) dr = 2F(x, 0)0(x, V). (4.4) 
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Equation (4.3) can therefore be expressed as 


KIN ( 


(P(x, t) ~ V,)+ (x, (4.5) 


2, 
where, from (4.4), V, and V, are given by 
—%.V, =n, +m, = (ge t+ | 
> (4.6) 
-%.V, =n, = x.U,—(gi + 
so that 
x.(U,+V;) = —Kxe(x), x.(U.+V,) = (4.7) 
c(x) being, as before, the speed of free surface waves of wave-number x. 
Notice that the velocities V (relative to the convected frame) of the frames 
of reference in which the time scales are to be measured are each specified 
by only one scalar equation, so that either the speed V of this frame, or its 
direction of motion can be chosen arbitrarily. For convenience the 
direction of V can be chosen to coincide with that of U, and so 


(U,+V,)cosz = —c(k), (U.+V,)cosa = c(«). 


Hence V, = V, = —U.4+¢(«)/cos«. (4.8) 
The absolute speeds of the frames of reference in which the scales @ are 
to be measured are c(«)/cosx and —c(«)/cos%, the negative sign implying 
motion in the direction against the wind. It is evident that the time scale 
of the pressure fluctuations observed in a reference frame moving against 
the wind is less than that moving with the wind at the same speed, so that 
A(x, 1,) < A(x, V,), and the expression (4.5) can be approximated by 


K2] I(x)t A(x, V). (4.9) 


2 2p?n2 

The speed I” = c(«) cos x is just equal to that of the waves of wave-number « : 
propagating in the direction x, observed along a line parallel to the wind, 
so that the integral time scale represents the time scale of the variations in 
phase difference between the pressure fluctuations and the wave pattern. 

‘The expression (4.9) is the basic result of this part of the paper, and 
from it several important deductions can be made immediately. Since 
the only way in which the elapsed time ¢ enters the right-hand side is in 
the explicit linear factor, the mean-square wave height is directly propor- 
tional to time in the principal stage of development. Furthermore, the 
shape of the wave spectrum is independent of time, at any rate until the 
development is such that the mean-square slope (governing the importance 
of non-linearities in the equation of motion) associated with any group of 
wave-numbers becomes too large. It will be seen later that this occurs 
first at the shorter wavelengths excited, and the contribution from these 
to the total mean-square surface displacement is small. Thus the direct 
proportionality of & to ¢ should continue until the mean-square slope of 
the largest waves attains a limiting value, after which time an equilibrium 


(D(x, t) ~ 


structure is presumably attained. 
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However, a direct application of (4.9) is hindered by the presence of 
the factor 6(x, V) about which little is known. When V = U,, so that the 
wave pattern moves with the convection velocity of the pressure fluctuations, 
we have the true resonance waves which figured in $3.2, and the time scale 
6(x, V’) is equal to the convected or natural time scale ©. The use of the 
capital letter in this latter case emphasizes the particular significance of 
this time scale. However, when V = U,, and usually V < U,, the time 
scale in the reference frame moving with speed I” is less than ©, since it 
is governed by the rate at which the convected fluctuations are swept past 
the point of observation. A rough approximation to 4(x, 1’), but one which 
is probably satisfactory, is to take it as the length scale «! divided by the 
relative speed of the convected pressure fluctuations and the waves; thus 


x)’ 


A(x, V’) (4.10) 
when U, is greater by about a factor of about two than c(«)/cosz. More 
properly, the time scale represented by the expression on the right is a 
differential scale, whereas 6(x, V’) is defined as an integral scale. However, 
the two scales should be proportional, and probably of the same order of 
magnitude, so that the approximation (4.10) can be used in the absence 
of better information, provided the results are interpreted as being valid 
only to orders of magnitude. 
For gravity waves, n3 = gx, so that, from (4.9) and (4.10), we have 


‘ (x)t 
06) ~ — ele) 


(4.11) 


When U is greater than c(«)/cos% by a factor of about three or four, this 
expression can be further simplified to 


I] (x)t 
D(x, t) ~ (4.12) 
so that we have immediately 
| 
Bw | 4.13 
|’ (4.15) 
and, for the mean-square slope, 
(ver ~ —CP) (4.14) 


20 2p?U.g 

Further reference will be made presently to these particularly simple 
and important results. At this stage we must consider the conditions 
under which the assumption that U, > ¢(«) cos is likely to be valid for 
most of the range of wave-numbers over which the spectral function (x) 
is significant. Support is given to the approximation (4.10) by the 
observation that for moderate winds over all but the longest durations, 
the wave spectrum is appreciable only over wave-numbers for ‘which the 
condition U,, > c(«) is satisfied. It is instructive, however, to see why this 
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is so by examining the relevant meteorological and other evidence on the 
surface pressure fluctuations. ‘The dependence of this function on x is 
determined, as (4.9) shows, jointly by 6(x, 1’) and the pressure term, which 
for gravity waves where 5 = gx is of the form «II(x). It is unfortunate 
that the properties of the spectral function I(x) of the imposed surface 
pressure distribution are at present almost unknown, though by using the 
results of some recent atmospheric measurements of the turbulent velocities 
(made over land, alas, not over the sea), some of its qualitative properties 
can be surmised. But one important property that can be established without 
recourse to these experiments is that when x = 0, I(x) = 0; this is proved 
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in the Appendix to this paper. ‘Thus in any direction « in the «-plane, 
Il(x) rises from zero at the origin to a maximum and then decreases again 
to zero as « approaches its value for the dissipating eddies. We are now 
faced with the question of deciding the order of magnitude of the values 
of « for which II(x) is a maximum. ‘The only guide we have at present 
are experiments such as those of MacCready (1953 a, b), Taylor (1955) 
and Ellison (1956) on velocity autocorrelations over land in winds within 
the range 1 to 5msec"!. ‘These point to a longitudinal scale of the energy- 
containing eddies in the lowest few metres of the order 10 m (the lateral 
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scales are considerably smaller), that is to say, the peak value of the velocity 
spectrum under these conditions occurs at wave-numbers of order 
27/108 =6x10-%cm™. Apparently this scale does not depend greatly 
upon wind velocity or distance from the ground. In isotropic turbulence, 
Batchelor (1951) has shown that the scale of the pressure fluctuations is 
rather less, about half, of that of the velocity fluctuations; and if a similar 
situation exists in the atmosphere, the maximum value of II(x) may occur 
for values of « of about 10-*cm~!. Furthermore, atmospheric turbulence 
generally contains a large range of eddy sizes, so that the spectra tend to 
be fairly broad. These considerations make it likely that I(x), in a given 
direction x in the «-plane, is a function of « of the shape illustrated in the 
top portion of figure +, with a broad maximum near« = 10-?cm', decreasing 
to zero for both large and small values of «. The function «II(x), being 
weighted towards the larger values of «, has its maximum at a larger «, and 
would be expected to be of the form also illustrated in figure 4. 

The time scale 6(x) has its maximum when « has such a value that the 
waves of this wave-number move with the convection velocity U.; then 
&(x) = ©, the development time of the pressure fluctuations. For a wind 
speed of, say, 10? cmsec™, a fresh sailing breeze, the maximum for angles 
in the range — 47 to 37 (containing most of the wave energy) occurs when 
k~ 10% cm. For larger and smaller values of « away from this maximum, 
6(%) is given by (4.10), so that the function has the form of the solid curve 
in the upper part of figure 4. 

The wave spectrum is given by the product of the functions #(«) and 
«I1(x) represented by the curves in this diagram, and it is clear that for a 
wind speed of 10msec~! the range of wave-numbers PP’, over which the 
approximation (4.10) is insufficient, contributes little to the integral { dx. 
The result (4.12), therefore, would be expected to give a good approximation 
to the wave spectrum except for the smallest values of «, whose contribution 
to the mean-square wave height is small. For wind speeds greater than 
the value 10°cmsec"! chosen for purposes of illustration in figure +, the 
wave-number for which U,. = c(x) decreases as U~*, so that the curve (x) 
is translated to the left. It is unlikely that the reciprocal of the scale of the 
pressure-generating eddies should decrease with wind speed as rapidly 
as this, so that the range of wave-numbers PP’ contributes even less to the 
integrated wave spectrum, and the approximations inherent in (4.12) and 
(4.13) become more accurate. However, when the wind speed and associated 
convection velocity is small, say 1 msec "|, it is probable that the range PP’ 
will occur at a point where the value of «II(x) is no longer negligible. “Vhese 
circumstances are illustrated in figure 5, and the wave spectrum calculated 
from the curves 6(x) and «II(x) is clearly no longer similar in shape to II(x) 
but possesses a more pronounced peak. In light winds, therefore, the 
wave spectrum is narrower, and presumably the wave pattern more regular, 
a deduction which would appear to be substantiated by visual observations. 
A similar situation exists for higher wind velocities when x is very nearly 
+47, and the curve 4(x) is displaced considerably to the right, but the 
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range of angles over which the approximation (4.10) fails is small and again 
contributes little to the integrated wave spectrum, or to €?. 


4.2. Termination of the principal stage of development 

There is another effect which limits the validity of (4.12) and (4.13) 
and which we have not yet considered. As the wave pattern develops, the 
mean-square slope of the waves increases linearly with time. After a certain 
interval, the mean-square slope associated with some band of wave-numbers 
becomes so large that the waves can no longer be described by a linear 
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theory. The non-linear processes which then become important tend, 
among other things, to limit the maximum surface slope that can be 
attained. Without a detailed examination of the role of the non-linearities 
in modifying the wave spectrum, it is not immediately clear exactly what 
criterion determines the importance of the non-linearities, although it is 
very reasonable to suppose that it concerns the value of the slope spectrum 
«*I1(%) in some manner. There is some evidence (Danel 1956) that the 


effect is appreciable when the mean-square surface slope (VE)? = [«x?@(x) dx 
attains a certain value. But whatever the exact criterion might be, provided 
only that it involves a limitation on the slope spectrum, which is largest 
when the wave-number « is large or the wavelength A small, its effect will 
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clearly be felt first in this range of «; the actual wave spectrum will for these 
wave-numbers be smaller than that predicted by the linear theory, and 
instead of continuing to grow will presumably attain some sort of equilibrium. 
As time goes on and the wave amplitudes at other wave-numbers increase, 
the wave-numbers contributing significantly to the slope spectrum extend 
to smaller and smaller values ; so that the effects of the non-linearity gradually 
spread to the wave-numbers for which the wave spectrum ®(x) itself is 
appreciable. ‘This process whereby the principal stage of development of 
each wave-number is ended in turn implies that the value of A for which 
(x) is a maximum increases with time, since the longer waves continue to 
grow while the amplitude of the shorter waves is limited. A detailed 
analytical description of this process has not yet been attempted, but it 
seems likely to provide part of the explanation of the observed increase 
with time in the dominant wavelength of the surface displacement. 

An important consequence of these remarks is that the result (4.13) 
remains valid until the wave-numbers near the maximum of the linear 
spectrum (x) become affected by the non-linearity, although over a large 
range of larger wave-numbers, or smaller wavelengths, the principal stage 
of development has terminated, and the expression (4.14) for the mean- 
square slope has become seriously inadequate. The reason is, of course, 
that the wave-numbers contributing greatly to the slope spectrum make 
little contribution to the wave spectrum. 


4.3. Comparison with observations on the wave field 

The most important result of §4.2 from the practical point of view 
is (4.13). In order to compare this result with the oceanographic measure- 
ments that are available for conditions of finite duration, we must obtain 
an estimate for p? in terms of the quantities that determine the mean velocity 
profile of the wind. It has already been pointed out that these relevant 
parameters are the friction velocity uw» and (for aerodynamically rough 
flow) the length z), so that the mean-square pressure fluctuations at the 
surface should be expressed as 


= Apiut 


(4.15) 
where A is a numerical constant to be estimated and p, is the density of 
the air. 

In shear-flow turbulence in general, p? is of the order of }p2 v?V?, where 
v and V are turbulent and mean velocities representative of the region 
concerned. ‘The atmospheric measurements of MacCready (1953b) are 
helpful in applying this expression and obtaining an estimate for A. Under 
conditions of almost neutral thermal stability over land with 2, of the order 
0-1cm, he found that in the lowest few metres v? = v7 +73+72 = JV? 
In view of the fact that under these conditions, the maximum of the pressure 
spectrum occurs at wave-numbers corresponding to length scales of order 
10 m, it seems likely that the velocity fluctuations at heights of the order 
of 1m contribute most to p*, and by relating MacCready’s measured 


q 
j 
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values of v*, V and uy, we obtain the result that, to the correct order of 


magnitude at least, 
= 9x (4.16) 


In most published oceanographic measurements, the wave heights 
observed have been expressed in terms of wind velocities U measured at 
some level, usually several metres, from the surface. In order to compare 
our theory with these observations, it is necessary to relate the friction 
velocity to the measured ‘wind velocity’ for conditions similar to those 
under which the observations were made. For the sea surface with moderate 
winds, a typical value of 3, is of the order of 0-15 cm, and if the height z 
at which U is observed is of the order of 5 m, we have from (2.1) that 
approximately 

U = 18uy. 
‘lhe factor 18 is not critically dependent upon z or 3,, so that it can reasonably 
be taken as representative of the conditions under which measurements, 
such as those reported by Sverdrup & Munk (1947), were made. From 
(4.16), therefore 
= 9x U*. 

This relation can be substituted into (4.13), and since for pressure components 
of the scale of gravity waves the convection velocity U, is equal to the 
velocity several metres from the surface, it is approximately equal to the 
‘wind velocity’ U. We therefore have, from (4.13), 


2U% 
0-035(£") 
Pw) 8 


where p,, is the density of the water. 

In oceanographic practice, wave heights are usually measured from 
trough to crest, and a widely used measure is H, the mean height of the 
third highest waves. Longuet-Higgins (1952) has shown that H? = 8 & 
approximately, so that the previous equation can be written 


o 1/2 


taking the value 10°° for p,/p,,.. The accuracy of the constant in this 
expression is quite low, mainly as a result of the uncertainties in the value 
of p® to be expected upon the surface. (Some reliable measurements of 
this would be very valuable.) It might be pointed out, however, that the 
factor 2-°* in (4.13) represents an underestimate in consequence of the 
neglect of 6(x,V,) and of the term c(x)/cosx in (4.11) leading to (4.12). 

When the result (4.17) is compared with some of the wave heights 
observed under conditions of finite duration as in figure 6, we find quite 
good agreement over the rather restricted range of durations for which 
measurements have been made. ‘The observations are those given by 
Sverdrup & Munk in their figure 7. The ratios of the speed C of the 
highest waves to the wind speed U are also plotted in figure 6, and in these 
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observations the ratio C/U is of the order of unity. Under these circum- 
stances, other wave-generating processes such as sheltering and the effects 
of variations in shear stresses may have become important, and, in any case, 
this theory would be expected to give an appreciable underestimate of the 
wave heights for the reasons mentioned above. In spite of this, it seems 
that the approximations used are sufficiently accurate to account for the 
order of magnitude of the wave heights observed. 
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Figure 6. 


We are now in a position to see rather more clearly the probable reason 
for the failure of Eckart’s theory to predict the magnitude of the wave 
heights generated by the wind. His less precise specification of the pressure 
distribution has ‘smoothed off’ the resonance peak of the response of the 
water surface, and it is the wave-numbers near this peak that can contribute 
largely to the wave spectrum at large durations. 


APPENDIX 
We wish to show that II(x) = 0 when x = 0, or equivalently that 


[ +e) dr = 0, 


the integration being over the whole surface with x fixed. ‘This result is 
established by joint consideration of the surface boundary condition (2.12) 
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relating the surface displacements to the surface pressure distribution, and 
the dynamics of the mass of air above. If (2.12) is multiplied in turn by 
dA"(x, t) and da(x, t), and the probability average taken to find the spectral 
functions in the manner of (2.3), we can put x = 0 ( and so n, = n, = 0) 
in the resulting equations. The result of these operations can be expressed 
as 


| &(x)E(x +r) dr = 0, 


| p(x)é(x+r) dr = 0, 
or if vs is the velocity of the fluid normal to the surface &, say, 


| v,(x)v,(x+r) dX(r) = 0, 
(A.1) 


| 
| p(x)v;(x +r) dX(r) = 0. | 


We now turn to the dynamics of the fluid above the liquid surface. It 
has been proved formally by the author in another context (Phillips 1956) 
that the integrated surface pressure covariances are related to the acceleration 
covariances, integrated throughout the fluid, by an expression of the form 


| = p, | de’ (A.2) 


where r’ = 7}, 73, 75 is a point in the interior of the fluid, and the integration 
on the right is throughout the volume occupied by the fluid. g(x) is a 
fluctuating quantity of zero mean whose precise nature is, in this expression, 
immaterial; it is introduced simply to provide integrable covariances with 
p(x+r) and 7,(x+r). Taking g(x) = p(x), the surface pressure at the 
point x, it can be shown to be a consequence of continuity that 


| | dr dr; 


is independent of But when = 0, ie. at the surface ° 
vanishes according to (A.1) so that the integral 


| | | +4’) dr, dr’, 


vanishes throughout any layer of large finite thickness above the surface. 
Hence, from (A.2), 


| = 0, (A.3) 


which is the result that we set out to establish. 


| 
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A cross-flow theory for the normal force on inclined 
bodies of revolution of large thickness ratio 
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In the Munk~Jones cross-flow theory for slender bodies of revolution 
(Munk 1924; Jones 1946), the cross force on an inclined body is obtained 
by replacing the three-dimensional flow by a non-steady two-dimensional 
flow, and by equating the cross-force to the rate of change of cross-flow 
momentum on a transverse lamina moving past the body with the free 
stream velocity Uy. The result obtained for the lift force LZ on an element 
of the body is, for small angles of attack «, 


dL/dx = \pU;(dA/dx) 2x, (1) 
where A is the cross-sectional area of the body, and, by integration 
L = ipU5 Ayia, (1a) 


where A, is the base area of the body. 

Although the fact is not generally recognized, this result is in poor 
agreement with exact calculations and with experimental data for finite 
bodies of revolution, and in order to develop a cross-flow theory for bodies 
of large thickness ratio, it is necessary to examine the assumptions of the 
simple theory. 

‘The main assumption is that the flow in the cross-plane is approximately 
two-dimensional. ‘Thus, in the three-dimensional incompressible potential 
equation (using wind axes) 

Perr + yy 0, (2) 
the first term must be small compared with the other two. By writing 
x =1X, y=sY, x =sZ, where s is a typical width and / a typical length 
of the body, and ¢(x, y, z) = Uy 1M(X, Y, Z), we may re-write (2) as 
2 
+Oyy = 0. (3) 
It is evident that the condition that the first term is small is that s?//? < 1. 
For compressible flows satisfying the Prandtl-Glauert differential equation 
for the perturbation potential : 


+ Pee = (+) 
where 8 = (1— M?)!?, the compressible flow is related to the incompressible 
flow past a body of revolution elongated in the x-direction by a factor £, 
and the condition for two-dimensional flow in the cross-plane becomes, 
for the compressible flow, f?s?//? < 1. For moderate values of 8, this 
condition will hold for surprisingly large vaiues of the body thickness 
ratio. (For a 20° cone, for instance, s?//? = 0-1.) 
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In order to calculate the cross-flow, it is necessary to solve the two- 
dimensional equation 
$yy +. = 90 (5) 
with ¢,, specified on the boundary. In the Munk—Jones approximation, 
this boundary condition becomes (using cylindrical coordinates (x, r,«) 
with the body shape defined by R(x)) 


d 
¢,(x, R,w) = U,cosx R(x) + Upsin «cos a, (6) 
or, by superposition for the cross-flow, writing ¢, = $,,+ ?2,, 
R,w) = Uyxcosw (7) 


for small angles of attack. 

This approximate tangency condition amounts to an assumption that 
the longitudinal component of the velocity on the body, which is the velocity 
with which the transverse lamina moves along the body, is approximately 
equal to the main stream velocity Uy. ‘This is manifestly untrue for thick 
bodies of revolution. For a 20° cone with a free stream Mach number 
of 1-3, for instance, the surface velocity is 0-8U,), and the longitudinal 
component is 0-75U). ‘This suggests that an improvement to the cross-flow 
theory may be made by employing a more exact tangency condition in which 
the main stream velocity is replaced by the longitudinal component of the 
surface velocity U(x). ‘The lamina then moves along the body with a 
velocity U,cos #, where @ is the inclination of the surface to the body axis. 

Thus (7) becomes 


R,w) = U,cosw«. (8) 
The lift force on an element of the body is then 
dL/dx = \pU%cos?6(dA/dx) 2a. (9) 


In order to illustrate the use of this result, a comparison is made in 
figure 1 with the exact calculation by Kopal (1947) of the supersonic flow 
past yawed cones of semi-angle 9%). Since, for such flows, the surface 
velocity U, is constant, equation (9) may be integrated to yield 

of U, cos 
A», 2a. (10) 
Uy 
The lift coefficient slope is then 
[Cu] L U.cos Iy\? 
= 0" Agar Uy 
where U, is obtained from Kopal’s calculations. ‘To convert the quantities 
tabulated by Kopal to the notation used here, the following relationships 
were employed: 
U, = U,{v? + 2a?/(y—1)}!2, and C1, = Cy,—Cp, 


where Cy, = 8K, 7, 8K 
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The results given in figure 1 for cones of semi-angle 10°, 20° and 30 
show that the present method gives good agreement with the exact results 
for cones up to 30° semi-angle, and correctly predicts the variation of Cz, 
with Mach number for quite large Mach numbers. ‘This implies that the 
cross-flow close to the body remains approximately two-dimensional even 
for cases in which the parameter {?s?//? becomes large. The Munk—Jones 
value Cx, = 2(radian)~! is also shown on the figure. 
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Figure 1. Variation of lift force slope with Mach number for supersonic flow past 
inclined cones. 


The method has also been employed by the author (Cox & Maccoll 
1956) for calculating the cross-force on yawed cones in fully developed 
cavitating flow. 
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The laminar and turbulent mixing of jets of 
compressible fluid. Part I Flow far from the orifice 
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(Received 18 Fanuary 1957) 


SUMMARY 


lor flows of jet type, the assumption of a coefficient of eddy 
kinematic viscosity in turbulent flow leads to the possibility of 
combining in one the equations for laminar and turbulent motion. 
An approximation to the solution of these equations is found for 
the flow of compressible fluid issuing from a narrow slit, far from 
the slit. ‘lhe stream function is expanded in a power series in 
squares of the Mach number. Bickley’s solution (1937) for the 
corresponding problem in incompressible flow is used to start the 
iterative process by which successive terms of the power series are 
obtained. In order to find an analytical form for the second term 
of the series, it has been assumed that the Prandtl number is 
unity, that the viscosity varies as the mth power of the absolute 
temperature, and that the stagnation temperature of the jet is the 
same as that of the surrounding gas. ‘The solution found differs 
only slightly from that of Howarth (1948) and Illingworth (1949) 
when laminar flow is considered; only the ‘change of scale’ 
effect (arising from a distortion of the coordinates in Bickley’s 
solution) is of importance. In turbulent flow the effect of the 
second term of the series is as important as the ‘change of scale’ 
effect. ‘The effect of compressibility on the width of the mixing 
region is discussed for both laminar and turbulent jet flow far 
from the orifice. 


INTRODUCTION 


The problem of the mixing of streams has received a good deal of 
attention, both theoretically and experimentally, in recent years. The 
interest has been in the extension to compressible fluid flow of the problems 
investigated for incompressible fluid during the years before the war. For 
jet flows there are two main simplifications of the real problem; one is to 
suppose that the jet issues from a narrow slit, in which case the field at 
some distance from the orifice may be examined on the basis of a similarity 
of velocity profiles in sections normal to the axis of the jet, and the other is 
to treat only the so-called ‘half-jet’ mixing, being the case when the two 
streams in relative motion are both semi-infinite laterally. Each of these 
problems is itself divisible into two further problems, according as the flow 
is supposed to be laminar or turbulent. 

In this paper it is first shown that the assumption of the existence of an 
eddy kinematic viscosity coefficient in turbulent flow, which is known to 


F.M. 2H 
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be reasonably satisfactory for jet flows, leads to the possibility of combining 
in one the equations for laminar and turbulent motion. Next, an approxi- 
mation to the solution of these equations is found for the flow of compressible 
fluid issuing from a narrow slit, far from the slit. ‘The method used is that 
developed by Pack (1954) for the study of axially symmetric jets of com- 
pressible fluids far from the orifice. ‘The stream function is expanded in 
a power series in squares of the Mach number. ‘The solution of Bickley 
(1937), for the corresponding problem in incompressible flow, is used to 
start the iterative process by which successive terms of the power series 
are obtained. In order to obtain the second term of the series in an analytical 
form, it has been assumed that the Prandtl number is unity, that the viscosity 
varies as the mth power of the absolute temperature, and that the stagnation 
temperature of the jet is the same as that of the surrounding gas. ‘lhe 
solutions for two-dimensional flow out of a slit obtained by Howarth (1948) 
and Illingworth (1949) are seen as particular cases of the solutions found 
here. ‘The problem of the ‘half-jet’ mixing will be considered in a later 
paper. 
EQUATIONS OF MOTION 


Let u, v be the velocity components (or mean velocity components in 
the case of turbulent flow) parallel either to Cartesian axes (*, Vv) in two- 
dimensional flow, or to cylindrical coordinates (x, y) with ~ measured along 
the axis of symmetry in axially symmetrical flow. Let the origin of 
coordinates be taken at the point at which the mixing begins. Let p be the 
density of the gas in the jet and 7 its absolute temperature. 

When the Reynolds number of the flow is large, the approximate equation 
of motion for jet mixing is the same as for boundary-layer flow at constant 


pressure : 
ou ou 1 oO Ou (1) 
ax ay) yr Oy \P oy)’ 


where 6 = (0) for two-dimensional flow and 6 = 1 for flow with axial symmetry. 
In this equation « = p/p, the ordinary kinematic coefficient of viscosity, 
when the motion is laminar, while when the flow is turbulent ¢ = e(x), 
a coefficient of eddy kinematic viscosity which is assumed to be independent 
of the y-coordinate. 
The equation of continuity is 
(puy®) + < (pey*) = 0. (2) 


From this equation it follows that a stream function y exists such that 


pu cy pr. CX 


If the independent variables are changed to (x, 2), where 


pratt? = (148) py*ay, (3) 
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the equation of motion becomes 
p, 2 dz \2° dz 


Oz dxdz® ds oz Os 
To allow for the variation of viscosity with temperature, write 


lad = ‘“*\n 
= (7) = (rey. (5) 

The boundary-layer assumptions underlying the equations of motion 
imply that pressure is constant, bay at least a first approximation, through the 
field, giving p*7* = 1, where p* = = p/p. 

Now write € = e(x), w ted (i) in a laminar jet, 6 = w/p and e(x) = 1, 
and (ii) in a turbulent jet, 9 = e9, Reichardt’s constant exchange coefficient, 
with e(v) an experimentally determined function of x. Introduce also the 
new independent variable 


Then the equation of motion reduces to 


where % = p, and 8 = n—1 for laminar flow, and « = ep, and B = —2 


for turbulent flow. 
‘The energy equation is 


1 0 (y*pep ou\2 
of Ox =) ( C p =) (; ~) (7) 


where €7, = «/p for laminar flow, « being the coefficient of heat conduction, 
and when the flow is turbulent €7 is the eddy coefficient of heat transfer, 
and where 7 = C,, T (the enthalpy per unit mass), C,, being the (supposed 
constant) specific heat at constant pressure. 

When it is assumed that the Prandtl number o, which is equal to C,, €/e,, 
has the value unity, a particular integral satisfying equations (1), (2) and (7) 
is the Crocco relation 


A+ Bu, (8) 
where A and B are constants, the values of which depend upon the boundary 
conditions in the problem considered. 


‘TWO-DIMENSIONAL MOTION 


For two-dimensional flow, 6 = 0 in (6). In the first place it is of interest 
to note the simplifications of the equation of motion which are then possible. 

When 7* = 1, it follows that z = y, p = p,; and the equation is that of 
incompressible flow in a boundary layer. ‘The same equation results when 
8 = 0, this corresponding to laminar flow of a compressible fluid with 
viscosity proportional to temperature (m = 1); for a two-dimensional jet, 
this case was treated by Howarth (1948) and Illingworth (1949). 
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FLOW FAR FROM THE ORIFICE IN A TWO-DIMENSIONAL JET OF COMPRESSIBLE 
FLUID 
The solution tor the flow far from the orifice in a two-dimensional 
laminar jet of incompressible fluid was given by Bickley (1937). It is 


“y= 3tanh a), (Y) 


where 7» = = @® and a is a scale factor related to the momentum flux out 
of the orifice. It is not difficult to show that, if F is the momentum flux 
per unit time out of unit length of the slit, then a® = p, F 482°. 

‘he same solution may clearly be used for the turbulent incompressible 
jet with the proper interpretation of the symbols. When y = yx, € = 2°?, 
the kinematic eddy viscosity coefficient is proportional to x!*, and the 
solution is identical with that given by Gortler (1942); in this case the jet 
spreads linearly with x, in accordance with the experimental results of 
Férthmann (1934). Géortler’s results were in excellent agreement with 
experiment except near the edges of the jet. 

For the jet of compressible fluid, the boundary conditions far downstream 


may be expressed as follows: 


on y= 0, 0, and = = 0): 

5 
(10) 

on. u = 0, — = 0), andp =p 
where 7 and p, are respectively the undisturbed temperature and density 
in the outer medium into which the jet spreads. 

When it is assumed that the stagnation temperature of the jet and the 
temperature of the surrounding fluid are the same, Crocco’s relation (8) 
takes the form 

+l hy, (11) 
where 7, is the stagnation enthalpy in the jet. 

Phe solution (9) is the starting-point for the solution of the equations 
of compressible flow in jets. It is treated as a first approximation to the 
required solution. First it is seen that, from (9), 

barx 


p,u = sech®an. 


Insertion of this value into the Crocco relation (11) gives for the temperature 


y= 


sech'an, 


 Whene x'* and 7, is room temperature, 
T* ditfers from unity by less than 1”,, for jets of air with Mach numbers 
onthe axis V7, | 0-22. For higher Mach numbers it will not be sufficiently 


where k? = 18a'x? p? T,. 


accurate to use the solution for the incompressible jet flow to vield quantities 
in the compressible jet flow. In this case the technique used by Pack (1954) 
for axially symmetrical laminar jets may be employed. This requires the 


setting-up of an iterative process similar to that of Jansen (1913) and 
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Rayleigh (1916) by which the approximate solution at any stage is used to 
obtain a better approximation 
The stream function y is to be expanded in a series of the form 


This, substituted into equation (6), leads to the following equation tor F), 
on equating coefficients of k?: 

F, +2tanhéF] +8 sech*éF, + 4tanh sech*éF, = 28 sech*¢(7 tanh*é — 1). 
where the accents denote ditferentiations with respect to € = ay. With 
t=tanhé, F, =(1-—¢f*)G and (1-—¢*)* dG dt = H(t), the above equation 
simplifies to 

d*H dH 

( dt* dt 
A particular integral of this equation is the following polynomial of the 
sixth degree in ¢: 

H = 4 8(17 —72t*? + 45t* — 142°). 
The complementary function is a linear combination of Legendre functions: 
H = AP, (t)+ BO, (t), 

where P(t) = 1; 3(1-—2)) 
and QO, (t) = F(—a,, 1; 3(1+2)), 
with n, and n, the roots of the quadratic equation n(a+ 1) = 4, the notation 
being the usual one for hypergeometric functions. ‘The boundary conditions 


to be satistied are: 
ont = 0, F;(0) =0 and F;(0) is finite ; 


ont = 1, lim (1 = 0, i.e. lim = 0, | 
tl 


and lim {(1—1)*F;(t)—(1-t)Fi(t)} = 0. | 


t>1 


When the singularities of the Legendre functions near t = 1 are taken into 
account, it follows after some analysis that 4 = B = 0. Thus, finally, in the 
second approximation, when F, is calculated trom the value obtained for H 


«12 tanhé— sech?¢ (126 


+17tanhé- tanhie ) |. 
3 


A check on the error involved tn the use of the second approximation shows 
that it is likely to be less than 1°,, for W, < 0-52. In the neighbourhood 
of M, = 1 the maximum error would be about 2°, for laminar flow and 
14”,, for turblent flow. These estimates are based on the assumption 
that F, is of the same order of magnitude as FY on the axis of the jet. 

Thus the second approximation will express the etfect of compressibility 
with fair accuracy up to sonic speed on the axis. For values of WM, greater 
than unity the solution still formally applies, but the possibility of the 
occurrence of shock-waves may render it invalid. 


= | tanh ¢ 
aS 


; 
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The values of F,, Fj for various values of ¢ are given in table 1. 
The correspondence between y and z (or €), to the first order in k?, 
requires only the solution for incompressible flow. In fact, 


k* 
ay= E tanh (1 tanhte) |. 


This relation is naturally the same as that obtained by Illingworth for 
the two-dimensional laminar jet, his solution corresponding to the first 
term alone in the expansion of ys, for reasons already explained above. 


t E tanh~'t | 57(—2/B)F,(€) | 57(— 2/p)Fi(€) 
0 0 51-000 
0-1 0-100 5-012 47-852 
0-203 9-490) 38-777 
0-3 0-310 12-919 24-871 
0-4 0-424 14-794 7-888 
| 0-549 14°643 9-863 
0-6 0-693 12-019 | 25-646 
| ()-7 0-867 36°437 
| 0-8 1-099 2-454 39-027 
0-9 1-472 15-54 29-720 
| 1-0 * 36-000 0 


Table 1. 


The ditference between the velocity profite obtained from the first 
term in the expansion of is, when z is expressed in terms of the physical 
ordinate y, and the velocity profile in incompressible flow (which comes 
from the same term with z = y) will be called the ‘change of scale’ effect. 
‘The above expression for y in terms of € and ¢ shows that this effect tends 
to sharpen the jet, i.e. to decrease the width of the mixing region for both 
laminar and turbulent jets; e.g. when M/, = 1-0 the ‘change of scale’ 
effect diminishes the width of the mixing region by about 20°,,—the width 
of the jet being defined as 2¢?* a. ‘The numerical value of F{ is at most 
about 0-107 for a laminar jet. Since k?/¢?3 = 0-475 for a perfect gas 
with y = CC, = 1-4, the effect of the term F; on the velocity profile 
of a subsonic laminar jet is negligible. ‘This is not true for the turbulent 
jet because the value of Fy on the axis is —17/19. Since Fis negative, 
then changes sign and decreases to zero with increasing €, this term tends 
to make the profile of the compressible turbulent jet broader than an 
incompressible turbulent jet with the same axial velocity. For laminar 
and turbulent flow far from the orifice, the net effect of compressibility 
is respectively to decrease and increase the width of the mixing region. 

The analysis shows that the second approximation found in this paper 
will be useful for flows at any subsonic velocity. Higher approximations 


| | | 
| | 
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would require attention to the terms neglected in deriving the boundary- 
layer equations from the full Navier-Stokes equations. ‘The neglected 
terms are of a higher order than those retained in deriving the second 
approximation under the limitation imposed on the solution for the 
incompressible jet, namely p, M, 7? > a. 


MOTION WITH AXIAL SYMMETRY 
When 6 = 1 it is possible to carry out an iterative process just as above 
which obtains a solution for a jet with axial symmetry. ‘This was done by 
Pack (1954) for a laminar jet. The equation (6) shows that the same solution, 
correctly interpreted, may also be used for a turbulent jet. ‘The solution is 


9 


where k? = 26'/C,, pi, = and the relation between and the 
momentum flux F across a section of the jet is given by 5? = 3p, F/16za. 
The function F\(y) is defined by means of a new variable t = 1/(1 + }5*y"): 
F(t) = ——— (135 — 3088 — 84") — | n+ 


A table of values of F, and F/ when 2 = 0-76 is given in Pack’s paper (1954). 
‘The same general conclusions as were found for the two-dimensional 
laminar jet for the effects of ‘change of scale’ and compressibility apply 
equally to the laminar axially symmetric jet. It is found, for axially 
symmetric turbulent jets, that both of the above effects are of the same 
order of magnitude, but just as with the two-dimensional turbulent 
jet the sign of 7 -'F{(m) on the axis is negative; as 7 increases 7! F{(n) 
changes sign and then tends to zero as »— %. This behaviour of 
7 1 F{(y) tends to broaden the velocity profile of the jet, but the ‘change of 
scale’ is more significant and the net result is that the velocity profile of the 
jet is narrowed as the speed rises, in contrast with two-dimensional flow. 


One of us (L. J. C.) is indebted to the Sir James Caird ‘Trust for a 
‘Travelling Scholarship which enabled him to take part in this work. 
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Turbulent transport of heat and momentum from 
an infinite rough plane 
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SUMMARY 

In the first part of the paper the dimensional laws governing 
the processes of heat and momentum transport from an infinite 
rough plane are assembled and their consequences set out., In 
the second part, the detailed equations for the turbulent energy, 
the mean square temperature fluctuation and the covariance of 
temperature and vertical velocity are used, together with some 
speculative assumptions concerning the dissipative action of the 
turbulence, to derive a series of relations between the turbulent 
intensities and the Austausch coefficients. One of these relations 
indicates that the flux form of the Richardson number cannot 
exceed a critical value which is about 0-15. It follows that in 
highly stable conditions the buoyancy forces have little direct 
effect on the turbulent energy balance, their action being primarily 
to cause a reduction in the scale of the motion and some change 
in its structure. 


1. INTRODUCTION 

The theory of turbulent flow near a heated or cooled horizontal rough 
surface is notoriously difficult, but in the past few years some progress has 
been made. ‘This has largely been a result of an increased understanding 
of the role of the various physical quantities in the problem and the more 
precise formulation of a theoretical model, which have enabled certain 
relationships to be put forward on dimensional grounds. ‘These restrict 
the possible forms of certain functions (e.g. the velocity profiles) and show 
that some of the empirical formulae hitherto used to describe them 
(e.g. Deacon’s (1949) law) cannot be correct. It is the aim of the first 
part of this paper (§ 2 to § 5) to collect together these dimensional arguments 
and display the correct form of the physical laws. Nearly all of the steps 
in this demonstration have been described by earlier writers or are implicit 
in their work; but individual papers have presented only fragments of 
the argument and these have frequently been associated with assumptions 
which are very unlikely to be correct (e.g. that the Austausch coefficients 
for heat and momentum are equal) ‘Thus, since it is only after the 
dimensional arguments have been completed that it is profitable to examine 
the mechanism of the flow in detail, the author hopes that any readers who 
are familiar with these arguments will excuse their repetition. 
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The character of the second part of the paper (§6) is different. In it 
an attempt is made to learn more about the physical processes controlling 
the transfer of heat and momentum from the equations of motion. In 
order to do this, the detailed equations for the turbulent energy, the mean 
square density fluctuation and the covariance between density and vertical 
velocity are written down. ‘These alone are not sufficient to solve the 
problem, and to make further progress it is found necessary to resort to 
some speculation and guess-work ; thus the results that are obtained, although 
based on assumptions which the author finds plausible, cannot be regarded 
as having any established validity. ‘These assumptions concern the dissi- 
pative action of the turbulence on the turbulent energy, the density 
fluctuations and the covariance of the density and vertical velocity. It is 
suggested that the rates at which the turbulence would begin to destroy 
these quantities, if the simple processes continually creating them were to 
cease to operate, vary with the stability in sucha way that their ratios remain 
constant. It is shown on the basis of this assertion that the ratio of the 
Austausch coefficient for heat to that for momentum approaches zero when 
the flux form of the Richardson number approaches a critical value which 
is considerably less than unity. It follows that in very stable conditions 
it is not through the energy balance that the buoyancy forces exercise their 
great influence on the properties of the turbulence. ‘The intensity of the 
density fluctuations is shown to be an important quantity, and that further 
investigation of the way in which it is controlled is desirable. 


2. ‘THE THEORETICAL MODEL 

The boundary layers of meteorology commonly differ from those of 
aerodynamics in that their thickness is great and the region that is of primary 
interest and accesible for measurement is relatively small and close to the 
surface. ‘The actual thickness and possible growth of these layers is often 
assumed to be of little importance, and it is customary to treat ihem 
theoretically as if the properties near the surface were independent of the 
state of the flow at great heights. Moreover, the turbulent shear-stress 
and heat flux are considered to be independent of height. ‘This corresponds 
closely to the actual conditions near the ground so far as the shear-stress is 
concerned, but unfortunately the complex effects of radiation do in practice 
lead to some change with height in the turbulent heat flux. ‘These are of 
dominant importance in stable conditions when the turbulent Austausch 
coefficient is very small, and possibly also very near the surface. In addition, 
it is seldom that the ground or sea is sufficiently uni‘ rm horizontally for 
advection to be negligible and a steady and homogeneous state to be 
established. For the present, theory must ignore these complications and 
consider a steady state in which the turbulent fluxes are both independent 
of height. 

It is now possible to specify a theoretical model. We consider an infinite 
uniform rough plane (zs = 0) which supports a fluid of great depth. 
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A stress puX is applied to the plane tangentially causing it to move in the 
negative x-direction, and heat is supplied to (or extracted from) the plane 
at a constant rate. We take our frame of reference to be fixed relative to 
the plane and confine our attention to the steady state which the flow in the 
vicinity of the plane is supposed to attain ultimately. In this model the 
mean velocity of the fluid (relative to the plane) increases indefinitely with 
height in an unknown way; and because of this, the indirect and rather 
unnatural specification of the flow is needed. The unproved assumption 
that the flow near the surface attains a steady state is important as it implies 
that the flow there is not influenced by that at great heights. If it is false 
and the flow near the surface does depend on the thickness of the boundary 
layer, the great complexity of the meteorological problem would indicate 
little prospect of deriving any theory at present. 

We shall consider the fluid to be incompressible, but the known 
equivalence (subject to certain restrictions) between the motion of an 
incompressible fluid and that of a gas, enables the results to be applied to 
the atmosphere as required. | Moreover, we shall assume that the fractional 
changes in density are so small that their effect on the inertia of the fluid can 
be neglected, so that they need only be taken into account in the gravitational 
terms. In other words, the fluid is treated as having uniform density but 
variable weight, as is commonly done in problems of this type. 


3. "THE ROLE OF 3, 
In conditions of neutral stability* (that is, when there is no heat flux) 

it is known that the mean velocity U at height 3 is given by the equation 

dU tts 

dz Fe’ (1) 
where & is a constant with an observed value of about 0-4. The integral 
of (1) is 

kU 


Ux 


= log(s/39), (2) 


where 3, 1s a length characterizing the roughness of the surface; 3, enters (2) 
only through the boundary condition determining the constant of integration 
of (1), and so a change in 3, merely adds a constant velocity to the flow. 
It would seem reasonable to suppose that such a change would not affect 
the turbulence in any way (except at heights of the order of s,). Even 
in stable and unstable conditions, (1) still holds near the surface and we 
are able to state quite generally that the only effect of a change in the value 
of 3y ts to superpose a uniform translation on the whole flow without modifying 
its internal mechanism. 


* The word ‘ stability ’ is used in this paper to indicate the static stability of an 
inviscid fluid with the same vertical distribution of potential density as that in the 
flow considered; it does not refer to the growth of small disturbances in the moving 
fluid. 
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This statement is so obvious and elementary that at first sight it seems 
strange that it has received little attention from meteorologists. ‘lhe 
explanation of this probably lies in uncertainty concerning the assertion that 
the flow near the surface is independent of that at great heights. When the 
theory of the whole depth of the boundary layer (Ekman spiral) is considered, 
even in its simplest form (see Ellison 1956), it is found that zy is involved in 
determining its thickness. ‘Thus, unless the assertion is correct, 3) is 
likely to affect the whole structure of the flow. 


4. DIMENSIONAL ARGUMENTS 

‘The statement of the last section has some immediate and far-reaching 
consequences. It implies that the only variables deternuning the internal 
structure of the flow are gp’w/p (where p’ is the density fluctuation and p 
the mean density), v (the molecular kinematic viscosity), « (the thermometric 
conductivity), wg and z. At the very high Reynolds numbers of the atmo- 
sphere the molecular quantities may be omitted where they occur separately, 
but it is debatable whether or not the Prandtl number v/« ought to be 
considered in relation to heat transfer (Batchelor & ‘Townsend 1956). 
It is omitted here, but since it is a dimensionless constant the omission has 
little importance at this stage. _ Hence, in neutral conditions all dimensionless 
groups such as u?/uz, v?/uz, w?/uz, etc., (where 
u, v and w are the components of the velocity fluctuation and p is the pressure 
fluctuation) must be constant independent of height. Covariances 
such as u(x, Yo, must be functions of only (x, 

When there is a heat flux, an additional length L = puj/gp’w is available 
and all the dimensionless groups will depend on the dimensionless height 
variable, Z = s/L) (Obukhov 1946; Charnock 1956). ‘Thus, using the 
hypothesis that even in non-neutral conditions the effect of a change in 2% 
is merely to superpose a uniform translation on the whole flow, we have 

1 dU 

Uy d(log Z) 
For small values of Z the buoyancy forces are negligible and f’ must approach 
the constant value 1/k which it takes in neutral conditions. ‘lVhus, the 
integral of (3) is 


f (log Z|, sgn Z). (3) 


+log S,, = kf(log Z,,sgn Z), ©) 


where S,, is Businger’s (1955) stability index defined by 
S,, = = 22/2, (5) 


* These variables can, of course, be combined in various ways, e.g. to give 
(x, —X»)/(z; +22), ete. The choice of the most convenient form is a matter of some 


delicacy. 
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and f(log!Z|, sgn Z) tends to log Z for small Z. It will be seen from (4) 
that the velocity profile can take only two forms, one for stable and one for 
unstable condition; the logarithmic law is the limit of these for small 
heights. 

‘The determination of the function f(log Z , sgn Z) in (4) is, of course, 
a major task for theory, but not one that can yet be accomplished. However 
physical reasoning can take us a little farther. In 1920 Richardson showed 
that the quantity which is now known as the flux form of the Richardson 
number (Deacon 1955) (and denoted here by Rf) is equal to the ratio of the 
rate at which buoyancy forces extract energy from the turbulence to the rate 
at which it is supplied by the shear-stress. ‘This result is a direct 
consequence of the energy equation tor the turbulence, (13). Rf may be 
expressed in the following equivalent ways: 


Ku Kué dz gp'wi Ww 


Rf = Ri Ky(dU dz)? puzdU dz pu; 
where M~ 3(dU ds) ind Ky dp/dz 


Richardson’s arguments make it clear that Rf cannot exceed unity over any 
substantial depth of fluid, so it must havea maximum value (critical value) not 
greater than this. It will be suggested shortly that Rferit is in fact 
considerably less than one, but some deductions can be made from its mere 
existence. ARfand Ay, are functions of Z, and physical intuition urges that 
they are monotonic functions. — In that case, Rf must approach its critical 
value when Z becomes large and positive. We have in the limit: 


Rferit kiky Ky gpwKy pus. (7) 


Hence in very stable conditions Aj, is independent of height and determined 
by the heat and momentum fluxes. Linear velocity profiles such as are 
predicted by (7) have been observed in the atmosphere (e.g. by Rider & 
Robinson 1951), but the agreement is probably fortuitous since it is likely 
that in the conditions of the measurements the temperature profiles were 
controlled by radiation and that the turbulent heat flux was not independent 
of height. 

If Rr, as well as Rf, has a maximum value, (7) determines Ay in terms of 
it; but there is no reason to think that this is the case. It seems more likely 
(especially in view of (20)) that turbulence can be maintained at large values 
of Ri. = Proudman (1953), following earlier work of Sir Geottrey ‘Vaylor, 
has analysed measurements of currents and salinity in the Aattegat and 
found appreciable transport of momentum by turbulence associated with 
values of Ri up to 10. As Rr increases, the ratio Aj,/Ay, must decrease 
in proportion; so, in our model, Aj, decreases indefinitely with height and 
is not determined uniquely by the fluxes alone in the same way as Ay. 
In order to tind Ay, and hence the density profile, it is necessary to know 
how the limit is approached. ‘That is, one must know how Ay Ky, 


approaches zero as Rf approaches Rferit, and how Rf varies with Z. 
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At the other extreme, in very unstable conditions (large negative Z), 
a state of free convection may be assumed to exist, so that Ay is independent 
of uy. It follows that the quantity H* defined by, 

‘wp 
He. (8) 
is constant. Measurements made by workers in Australia (Taylor 1956) 
suggest that it has a value of about 0-8. An alternative expression for H*, 
which brings out the dependence of Ay on Z in these conditions, is 
7-1 ( 
H® = (9) 
Any sort of mixing length theory would seem to suggest that Ay, is also 
independent of wy 1n free convection, and this is generally assumed to be the 
case. It follows that Ag/Ay, must approach a constant value for large 
negative Z. 

Our knowledge of Ky as a function of Rf can now be summarized as 
follows. When Rf = 0, Ky, =k (see (1)); when Rf approaches Rferit at 
large positive Z, Ky, vanishes (see (7)); when Rf is large and negative, Ky 
is proportional to Rf!/4 (see (9) and (6)). Hence the well known and useful 
semi-empirical formula of Holzman (1943), 

(K y/k)? = 1—sRi, (10) 
where s is a constant, cannot be correct except for very small values of Rv. 
However, it is very convenient to have a semi-empirical formula like this, 
and so the following may be suggested : 

= 1— Rferit- (11) 
Note that the right hand side of(11) contains Rf in place of Ri. It can be 
written as a relation between Aj, and Z, since by (6) we have 
= 


(11) has no theoretical basis beyond having the correct form in limiting 
conditions and presumably a more elaborate formula with more empirical 
constants will be needed as soon as precise observations become available. 
In the case of free convection (11) implies 

KylKy = (12) 
When one inserts the numerical values already mentioned for A* and k 
and the value 0-15 for Rferit as is suggested below, this gives Ay/AKy = 1-6, 
which seems quite reasonable. 


5. "THE (LOG + LINEAR)-LAW 
Presumably A¥, can be expanded in a power series in 7: 


so RU /uy = log(z/2))—a, 7+... . (13) 
Businger (1955) has shown from an analysis of Rider’s (1954) velocity 
profiles that a, is about 0-8, while Monin (undated) gives a, = 0-6, 
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The semi-empirical equation (11) yields a, = k/4Rferit, which implies 
that Rforit = 0-14 if a, = 0-7. ‘This value has no more authority that (11) 
itself, but it does strongly indicate that Rfcrit is small. It also turns out 
that this value is virtually identical with that suggested by a quite separate 
argument in the next section. 


6. "THE TRANSPORT MECHANISM 

In this section an attempt is made to approach the subject from a new 
angle and to derive useful information from the equations of motions 
together with some relatively weak and general assumptions. ‘The 
suggestion already made that Rf has a small critical value is confirmed, 
and it is shown that the assumptions lead to a reasonably consistent physical 
picture of the processes governing the flow. It is not claimed that the 
validity of the assumptions is established; but, whether they are correct 
or not, the theory does indicate subjects which are likely to repay further 
investigation. 

From the Navier-Stokes equations, the equation of continuity and the 
equation of heat conduction, it is elementary to deduce the following 
equations for the mean square density fluctuation, the turbulent energy 
and the density flux: 


] Cp” ] dp 
le = — —duU 
55 q*) + + uw 
lcwp' wp'g —— 
WuV2u+vV20 + wV2w) = 0, (15) 
p 3 p 
rg ae — vp'V2w = 0 (16) 


In our model conditions are steady and so the first terms in these equations 
vanish; they have been written down merely to aid understanding. ‘The 
‘ditfusion’ terms like ¢(wp'*)/éz, etc., are probably always negligible since 
p, u2, etc., vary only slowly, if at all, with height. The‘ pressure-flow’ term 
0(pw),/cs in (15) is zero in neutral conditions, and there seems no reason to 
think that_it ever becomes significant. ‘There remain the dissipation terms, 
the term p’cp es in (16), and the simple terms of obvious interpretation among 
which we wish to discover new relations. 

Now it is well known that at high Reynolds numbers the rates of 
dissipation are very largely determined by typical length and velocity 
scales of the turbulence and not by molecular quantities such as viscosity. 
Let us therefore formally introduce decay times and 7, for 
u2 4024+ *, and wp’, such that, in the (imagined) absence of the producing 
terms, the mixing action of the turbulence would begin to destroy these 
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quantities at rates equal to 1/7), etc. In other words, let us define T,, T, 
and 7’, by the following equations : 


—dp 
wp 0, (17) 
—dU  — 
ut + weg = 0, (18) 
wp’ 
r, p (19) 


Since the mixing action of the turbulence tends to destroy the correlation 
between w and p’, as well as w* and p’, 7; may be expected to be substantially 
smaller than 7 and 7, which may well be roughly equal. While 7, 7, and 
T;, certainly vary with stability, we may reasonably guess that they are all 
dominantly determined by the same typical length and velocity scales and 
that the relation between the processes responsible for the destruction of 
the turbulent energy, the density fluctuations, and the density flux, change 
only little with stability. If so, 7), 7, and 7; must remain in roughly 
fixed ratios. Such an assertion is admittedly speculative, but seems to 
lead toa consistent physical theory. 

In (19), the term p’cp/ez has been incorporated in wp’/T;. This is 
plausible if the term is considered to represent the drag on a rising blob of 
hot fluid in unstable conditions; in stable conditions the term turns out to 
be negligible and so the interpretation is unimportant. 

It may be noted at once that (17) alone implies that Kg is always positive 
and that no selective action of the buoyancy forces, such as was once suggested 
by Priestley & Swinbank (1947), can cause a transfer of density against 
the gradient in stable conditions. 

After a little algebra, the equations (17) to (19) can be rearranged to give 
Ky/K yy T2 (say), and p” in terms of the known fluxes, the ratios of 7,, 7, 
and 7, and the cecpenants of the turbulent energy. Let us first examine 
the equation for K 

Ky + w*)w*{1 — Rf(1 + +2? + 7, 

Ky RIP 
which is the basic result of this part of the paper. There can be little doubt 
that for small positive values of Rf the term in in square brackets is the dominant 
one. If we take 7,/T, = 1 and (uw? +7? + w?)/w® to have the same value as 
in neutral conditions, which is known to be roughly 5-5, we see that Ky/Ky, 
falls to zero when Rf = 0-15. This small critical value of Rf is most 
remarkable since it implies that in stable conditions the buoyancy forces 
do not have any great effect on the energy balance. In these circumstances 
it is difficult to believe that the assumption that (uw? +7? + w?)/w? has the 
same value as in neutral conditions is far wrong; but, in any case, one would 
expect that in stable conditions horizontal motions would be favoured 
at the expense of vertical ones, thereby causing the quantity to increase 
rather than decrease. Such a change would lead to an even smaller value 


(20) 
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of Rfcrit. It is, of course, conceivable that 7/7, decreases in stable condi- 
tions, but there is no obvious reason why this should be expected, and it 
would require a very drastic change to affect the general conclusion that the 
critical value of Rf is small. 

‘The author knows of no measurements which are nearly precise enough 
to enable the conclusion to be checked. ‘The work of Proudman (1953), 
which was referred to earlier, indicates a value around 0-25; but the 
conditions of the measurements were very different from the ideal considered 
in the theory and the discrepancy is not disturbing. 

Even in neutral conditions it is known only that Ay/ Ky is roughly 
unity. ‘The work of Rider (1954) suggests a value about 1-3, while that of 
Swinbank (1955) suggests about 0-7. Cramer & Record (1953) found that 
K yy was often less than Kg, in agreement with Rider, but their results were 
detinitely affected by the trajectory of the air in which the measure- 
ments were made. If we take 7/7; =6, w?/uZ = 1-6, and again 
(u? +w?)/w? = 5-5, in neutral conditions (Rf=0), (20) gives 
Ky/Ky = 1:2. ‘That it should be necessary to take quite such a high value 
of 7T,/T; is perhaps a little surprising, but quite unavoidable. It seems 
likely that with these values for the constants, (20) will give a good representa- 
tion of Ky/Ky while Rf is small. 

With increasing instability the ratios of the intensities will begin to 
change. At the extreme, in free convection, the scale of the motion is 
limited to be proportional to the height; but the turbulent energy is produced 
entirely by the buoyancy forces, and so, being independent of the shear-stress, 
must on dimensional grounds increase as Z?%.  ‘Vhe equation for Ay/ Ky 
now degenerates into one for Ay alone, genie 


+72 + w*)w*{1 + (u2 + + w)/T, w?] 
Ky = 21 
H 2g" pw ( ) 
But by (9), 
= H®3 7 = = 
so 
(sgp'/p)'? = T, za (22) 


which connects the turbulent energy with the heat flux in free convection, 
We do not know the value of (w?+77+ a? 2)? and the ratio of the 7’s in 
these conditions. If (22) is considered as an equation for u? +7? +w?, 
it is not very sensitive to the former and we may take a value of 3 as being 
reasonable; in the absence of any evidence that the ratios of the 7”s differ 
from their values in neutral conditions, we can only make a provisional 
guess that they are the same. If so, we obtain 


+ w = 3(agp'w/p)?*. (23) 


Now let us turn to the second of the equations derived from (17) to (19), 
that for 7',, which is just a restatement of (18). 7’, has little direct interest 
since it is not measurable, but on account of our assertion that it is determined 
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by typical length and velocity scales of the turbulence, it can to a limited 
extent help us to form a physical picture. (u2 + 0? + w2)!2 2 is the obvious 
choice for the velocity scale, and so the equation for 7’, can be used to deter- 
mine a length Ly defined as 7,(u2 +7? +2? sla From (18) and (6) we 
obtain 

Ly = (v2 +7? + Rf). (24) 
This implies that 1, is proportional to z in neutral and unstable conditions 
and to L in very stable conditions, as was perhaps to be expected. ‘The 
approximate numerical values are as follows: when Rf = 0, Ly = 10z; 
in free convection, Ly = 5z; in very stable conditions, Ly =5L. It is 
of course only the relative sizes of Ly in the three states that is significant. 
It is a little surprising that its value in free convection is only half that in 
neutral conditions. ‘This result may be partly spurious and due to an 
incorrect numerical constant in (23), but it is hardly possible that that is 
the whole cause of the difference. It is more likely that the effect is in some 
way connected with the fact that in free convection the turbulent energy 
is more nearly isotropically distributed than in neutral conditions. 

The remaining equation to be discussed is that for p’?, 

(wp)? 1—Rf[1+ +02 +w*)/T, 

pw? 2(7,/T;)(1— Rf) 
If we keep the same numerical values as previously, this predicts that the 
correlation coefficient, wp’ (p'2)!3 2(w?)!?, is about 0-6 in free convection; 
0-3 in (almost) neutral conditions; and vanishes in very stable oiauiidiaias 
In principle these predictions can be tested easily, but again adequate 
measurements are lacking. ‘The values given by Swinbank (1955) are 
very scattered, but indicate a value of 0-4 rather than 0-3 in nearly neutral 
conditions. ‘The measurements do not show the form of the variation 
with stability beyond the sign of the general trend. 

The simple physical explanation of the relative inefficiency of heat 
transfer in stable conditions expressed by (20) and (25) is that then a 
displaced fluid particle tends to return to its equilibrium level before it 
has had time to mix with its surroundings. <A particle can transfer 
momentum during a brief excursion without mixing with its surroundings 
through the agency of pressure fluctuations, but in order to transfer heat it 
must mix. 

Some idea of the vertical distance travelled by particles before either 
returning towards their equilibrium level or mixing can be obtained from 
the length Ly = (p')!2/dp/dz. It can easily be shown that 

L2 2 1— Rf/Rferit 2% 
272 1—Rf* ’ ( ») 
which shows that in stable conditions Ly becomes much smaller than Ly,. 
No precise interpretation of this can be given since the two scales have a 
different nature, but it presumably means that in stable conditions scales 
measured in the vertical direction are liable to be much smaller than those 
measured horizontally. 


(25) 


F.M. 
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One last point is worth mentioning. In very stable conditions the last 
two terms of (19) become much greater than the first, and balance. Thus 


pdp/dz 


p dz 

In other words, the scale Ly is just such that the work required to lift 

a particle a distance L,, from its equilibrium level is equal to its average 

kinetic energy. 

7. ‘TRANSPORT OF INERT POLLUTANT 

If the theory of the previous section is applied in a similar way to the 
transport of an inert pollutant (i.e. one that does not affect the density 
of the fluid), it is found that the Austausch coefficient for it is necessarily 
equal to K,. ‘The reason for this is that no mechanism is included in the 
theory which would enable the turbulence to destroy the correlation between 
density fluctuation and pollutant concentration. ‘Thus the conclusion 
is one which may be seriously upset as the role of molecular processes in 
turbulent flow becomes better understood. Observational evidence based 
on direct measurements of the Austausch coefficients is at present conflicting ; 
and it is likely that measurements of the correlation between density 
fluctuation and pollutant concentration would provide a simpler and more 
sensitive test. 
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SUMMARY 

This paper deals with an experimental investigation of a 
turbulent low-speed jet of air spreading out radially over a flat 
smooth plate: a flow which has been discussed by Glauert (1956) 
in his theory of the wall jet. 

The aim of the experiments has been to determine the mean 
velocity distribution and rate of growth of the jet. It is found, 
within the experimental range and accuracy, that the velocity 
profiles are similar and that the rate of change of velocity and 
width of the jet can be expressed by simple power-laws. 


1. INTRODUCTION 


The term ‘wall jet’ was introduced by Glauert (1956) to describe the 
flow that develops when a jet, consisting of a fluid similar to that of its 
surroundings, impinges on a plane surface and spreads out over the surface. 
Glauert studied such a flow in two dimensions and in three dimensions 
with axial symmetry, and pointed out that it has features common to both 
the free jet and the ordinary boundary layer; thus, the spreading fluid is 
retarded by frictional resistance of the wall and the inner part of the flow 
may be expected to show a certain structural similarity to a boundary layer, 
whereas entrainment of quiet fluid occurs near the outer edge of the flow 
which accordingly is likely to resemble a free jet in character. 

‘This idea of a hybrid structure led Glauert to a solution of the turbulent 
wall jet by introducing an eddy viscosity distribution near the wall consistent 
with the Blasius power-law velocity-profile, and a constant eddy viscosity 
in the outer part of the flow. From this solution, which involved one 
disposable constant, it was concluded that exact similarity of the flow at 
all distances x from the origin—the ‘ point’ of impingement of the jet— 
could not exist, since the eddy viscosities in the inner and outer parts of the 
flow varied in a slightly different manner with x or, what amounts to the 
same thing, with Reynolds number. As a consequence, if the local ratcs 
of change with x of maximum velocity and jet width are expressed by simple 
power laws, the exponents of these laws will themselves vary slowly with x. 
However, it appears from the analysis that for any particular jet the variation 
is so slow that it would be undetectable experimentally. 

The experiments described in this paper were made with a turbulent, 
round, low-speed jet which spread radially over a smooth, plane surface. 
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Velocity distributions were measured and, as foreshadowed by the theory, 
were found to be similar, within the limits of experimental accuracy; they 
were also consistent with the shape of profile predicted by the theory. 


2. ARRANGEMENT OF THE EXPERIMENT 


The experiment arrangement is shown diagrammatically in figure 1. 
The wall was a flat plate of bakelite 125 em x 125 cm, placed horizontally 
on a steel frame. 


trom centrifugal blower 


| 28.4 
Range of experiments 


Figure 1. Arrangement of apparatus. All lengths in mm. 


The air was supplied by a centrifugal blower connected to a rubber 
tube of internal diameter 56 mm and length 2:8 m, which terminated, 
by way of a smooth contraction, in a pipe of 28-4 mm diameter and 1-07 m 
length. A flange of 126 mm diameter was attached to the end of the pipe. 
The reason for installing this flange was partly to prevent the generation 
of large scale disturbances at the pipe exit, and partly to get a closer 
approximation to the form of source assumed by Glauert. During all 
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experiments the pipe was placed at right angles to the bakelite plate with 
the exit 15-0 mm above the plate. Preliminary measurements showed a 
fully developed turbulent velocity profile at the pipe exit with a peak velocity 
of 34 m/s. With a previous arrangement in which the rubber tube had a 
90° bend the velocity profile at the pipe exit was not symmetrical. ‘The 
deviation was consistent with a secondary flow due to the bending of the 
tube and disappeared when the rubber tube was straightened. 

The velocity distribution along the plate was measured by means of 
a 1mm external diameter total-pressure tube connected to a ‘Casella’ 
U-tube manometer, the measuring accuracy being 0-01 mm water. The 
time required to obtain a reading with this accuracy was about 15 minutes. 


3. EXPERIMENTS 

The velocity distributions were measured at distances of 143 to 303 mm 
from the origin. At each station a Pitot traverse was made perpendicular 
to the plate and the velocities were computed by assuming the static pressure 
to be atmospheric everywhere. Each of the traverses consisted, on the 
average, of 20 readings at about 0-3 mm intervals near the wall and intervals 
of about 1-5 mm in the outer part of the flow. ‘The accuracy of setting the 
Pitot tube was 0-05 mm. 
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Figure 2. Experimental velocity profiles. 


In the range mentioned the maximum velocity varied from 6-19 m/s 
to 2-60 m/s. The width of the jet, measured from the plate to the point 
where the velocity had decreased to half the maximum velocity, varied 
from 11 to 23 mm. 

Figure 2 is a dimensionless plot of velocities obtained from nine traverses. 
The horizontal axis indicates the height vy above the plate measured in terms 
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of 6, where U(6)=3U,, and U,, is the maximum velocity. The vertical 
axis gives the measured velocities, U, in terms of U,,,.. Points corresponding 
to values of U/U,, below 0-2 in the outermost parts of the flow have been 
omitted. With the previously quoted estimate for the accuracy of a pressure 
measurement, the percentage velocity error is 7-5/U?, where U is measured 
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Figure 3. Variation of width of jet with distance from origin. 


in m/s. For the largest velocity presented in figure 2, 6-10 m/s, the error 
is 0-2°.,; for the smallest velocity, which occurs at the outer edge of the 
How at a point furthest from the origin, the accuracy of velocity measurement 
is 28°... Apart from errors in pressure measurement, there is a source of 
error at the outer edge where U tends to zero; the transverse inflow velocity 
tends to a finite value which results in a yawed flow relative to the Pitot tube. 
A rough estimate indicates that there is no yaw effect at UU), = 0-5, whereas 
at U/l’, =0-2 the true values are of the order 20°,, less than the plotted 
ones. ‘The points within a few Pitot-tube diameters from the wall are also 
subject to error. 

Figures 3 and 4 respectively are plots of log U,, (U,, in m/s) and log6 
(5 in mm) vs logx, where x is the distance from the origin in mm. The 
simple power laws quoted in the figures correspond to the least-square 
estimated inclination of a straight line through the points. ‘The estimated 
accuracy is +0-03 for —1-:12 and + 0-02 for 0-94. 
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Returning to figure 2, the experiments show no detectable departures 
from similarity of the velocity profile. Whether this observation holds for 
larger distances from the origin was not investigated due to the limitation 
set by the measuring accuracy of the apparatus. 
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Figure 4. Variation of maximum velocity of jet with distance from origin. 


4. COMPARISON WITH GLAUERT’S THEORY 

The shape of the velocity profile given by Glauert depends on a single 
parameter 2, which is uniquely related to «R'4. « is an empirical constant 
related to the eddy viscosity in the outer layer, and R (= U,, 6,/v) is the jet 
Reynolds number, where 6, consistent with Glauert’s definition, is the 
distance between the points at which U = U,, and U = }U,,. The relation 
between « and «R!! is given in Glauert’s table 1. 

Figure 5 shows the experimental profile and Glauert’s solution for 
(«R'4=0-102). This value of was chosen to give the bestagreement. 
with experiment. ‘The jet Reynolds number of the experiment was 3500. 
This corresponds to « = 0-013. ‘This figure differs slightly from the value 
quoted by Glauert, « = 0-012, which was based ona preliminary experiment 
at R = 5000. The deviation between the experimental velocity profile and 
Glauert’s solution is negligible except at the outer edge. However, this 
is analogous to the case of a free jet where the assumption of constant eddy 
viscosity leads to too large velocities in this region. ‘The experimental 
points closest to the wall are certainly in error, and should be discounted. 
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Corresponding to « = 1-3, the predictions given by Glauert’s table 1 
are that U, cx!" and 6 x x, whereas the present experiments give 
for these similarity exponents, — 1-12 + 0-03 for U,, and 0-94 + 0-02 for 6. 
The agreement in the case of U,, is reasonable. For 6 on the other hand, 
Glauert gives values greater than unity for all values of «. 
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Figure 5. Comparison of M. B. Glauert’s theoretical velocity profile and experiment. 


In the present experiment the Reynolds number varies as x°!’. By 
Glauert’s theory this would increase the value of y/6 at the velocity maximum 
by 0-003 for x varying from 143 to 303 mm. Considering the experimental 
accuracy, one could not expect to detect such a small change of velocity 
profile. 


The author is indebted to Mr M. B. Glauert, Mr H. K. Zienkiewicz 
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SUMMARY 

Detailed measurements have been made of the temperature 
field in natural convection above a heated horizontal surface in 
air, with and without a cold upper boundary. ‘The object was to 
provide experimental material for the testing of existing theories of 
convection and the development of new ones. Using a variety of 
experimental techniques, it was possible to measure the heat 
transfer, the mean temperature profiles, mean squares of the 
temperature fluctuations and the autocorrelation functions of the 
temperature fluctuations. ‘These results are considered in 
relation with the ‘similarity’ theory of Priestley and the ‘ neutral 
stability’ theery of Malkus. Both these theories lead to the 
conclusion that, for convection between parallel planes at high 
Rayleigh numbers, nearly the whole of the mean temperature 
variation occurs in comparatively thin surface layers whose 
mutual interaction is small. ‘The present experiments confirm 
this in some detail, but the exact form of the variation of mean 
temperature cannot be completely reconciled with either of the 
theories. 


1. INTRODUCTION 

The motion set up in a fluid which is enclosed by two plane horizontal 
boundaries and heated from below forms a very suitable subject for both 
theoretical and experimental study of heat convection, the boundary 
conditions being homogeneous in the horizontal directions. As _ the 
difference of temperature between the surfaces increases from zero, an 
ordered, cellular motion usually appears at a Rayleigh number about the 
predicted value of 1708, but this regular motion disappears and is replaced 
by a random ‘turbulent’ motion at some Rayleigh number around 50,000 
(Jakob 1949; de Graaf & van der Held 1953). More recent work by 
Malkus (1954 a) has shown that there may be as many as five transitions 
between one flow regime and another before fully turbulent convection is 
established for Rayleigh numbers exceeding one million. 

Apart from some observations of the motion of suspended particles 
(Malkus 1954 a), there have been few attempts to study any other aspects 
of the problem than the dependence of heat transfer coefficient on Rayleigh 
number, and the experimental work to be described was undertaken in the 
hope of providing more detailed information about conditions within the 
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fluid. Working with air, it has been possible to measure mean temperatures 
and temperature fluctuations at points between the surfaces and so to make 
a more stringent test of the several theoretical predictions than would be 
possible with heat transfer measurements alone. Similar measurements 
have been made in the convection field above a single heated plate. This 
last arrangement is an approximation to the convection over an infinite 
horizontal surface and may be compared with the predictions of the 
similarity theory (Priestley 1954). It may also resemble the convection 
near the lower of two parallel planes at very high Rayleigh numbers. 


2. "THE EQUATIONS OF MOTION AND OF HEAT 

Let us consider the convection in a perfect gas in the space above a 
rigid horizontal surface whose temperature is everywhere 7, measured in 
the absolute gas scale. ‘The coordinate system is chosen so that the origin 
is in the plane of the surface and Oz ts vertically upwards. Sometimes the 
Huid is bounded above by a similar horizontal surface in the plane z = D, 
and with temperature 7. If there is no such upper boundary, the tempera- 
ture at a great distance from the surface is 7. Then 


u,v, w are the velocity components parallel to Ox, 
Oy, Oz, g? = u? +07 
fT are the mean pressure and the fluctuation 
about the mean, 
T, 4 are the mean temperature and the fluctuation 
about the mean, 
y is the ratio of the specific heats, 
v is the kinematic viscosity, 
K is the thermometric conductivity, 
H is the constant upward flux of total heat, 
O = H(pc,, T) = (y~1)H yP is a ‘kinematic’ measure of the heat flux, 
A = gD log(T,/ T.)/v« is the Rayleigh number (this definition is 


consistent with use of the logarithm of the 
absolute temperature and, for numerical 
purposes, is nearly identical with the ordinary 
definition). 
For motions on the laboratory scale, and with fluid velocities small 
compared with the velocity of sound, the constant flux of total heat is given by 


Ary. 


pc, + pc (2.1) 


Cs 
where the single bar denotes the mean value with respect to time. This 


may be written as 


where QO 1s also a constant. 
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To the customary approximation, the equation for the turbulent energy 
per unit mass is 
éfi—-_—s dw 
w) = y 
(54 ) 87 (2.3) 


OZ \S 


where € is the rate of dissipation of energy by viscous stresses. ‘The terms 
on the left represent transport of energy from one part of the flow to another, 
the terms on the right represent generation of energy by buoyancy forces 
and the viscous dissipation. 

Another mean value equation of some interest is the equation for the 
entropy associated with the temperature fluctuations, which is — }c,, 6 /T?. 
To this approximation it is 


Ow 1 eT (Pw 
c FT T 3s (a3) Key (2.4) 
showing that entropy (or temperature fluctuations) is produced at a rate 
proportional to the mean temperature gradient. 

It may be noticed that all three mean value equations (2.2 to 2.4) involve 
temperature variations only as fractions of the local mean temperature, 
which suggests that they would be more naturally expressed in terms of the 
logarithm of the temperature rather than the temperature. ‘The distinction 
is of no importance if the total temperature variation is negligible compared 
with the absolute temperature, but this is far from true in most of the experi- 
ments to be described. In the following analysis of the experimental 
results, the logarithm of the temperature has been used consistently with 
a considerable improvement in the correlation with theoretical predictions. 


3. ‘THEORIES OF THE CONVECTION 

A first step in the study of any problem of fluid motion is to use 
dimensional analysis to obtain the forms of the functional relations between 
the observables. For any convection between parallel, horizontal surfaces, 
the boundary conditions are geometrically similar and may be specified by 
the separation of the surfaces and the difference of the logarithms of the 
temperatures. If the equations of motion, heat, and conservation of mass, 
are put into non-dimensional form by use of scales of length, temperature 
and velocity derived from D, log(7, 7.) and the quantities g, v, «, it Is 
immediately evident that the convection depends only on the Rayleigh 
number and the Prandtl number, vx. For example, the heat transfer 
coefhicient is 


OD 
and the distribution of mean temperature is . 
v 
— v 7 
log(T/7;) (5-42): (3.2) 


For very large separations of the surfaces (high Rayleigh numbers), 
it seems possible that the convection close to each surface may be nearly 
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independent of the presence of the other surface and so of D. ‘The basic 
equations may still be made non-dimensional, and conveniently so by using 
the following scales of length, velocity and (logarithmic) temperature, 


14 O02 1/4 
= (5) Uy = (ge)! 4, 


and their solution depends only on the Prandtl number. It follows that 


3 
log = (0, = (3.3) 
that log T. = *) (3.4) 
and that = (= (3.5) 


where C is a function of Prandtl number, and 7\, is the asymptotic 
temperature far from the lower surface but not so far as to be within the range 
of influence of the upper surface. If this temperature exists, equation 
(3.3) defines the variation of heat transfer coefficient with Rayleigh number 


cys Aus (3.6) 
log( 7/72) K 

for Rayleigh numbers so high that the temperature gradient is negligible 
except close to the surfaces. The proportionality between heat transfer 
coefficient and cube root of Rayleigh number is well confirmed by experi- 
ment for Rayleigh numbers above 108 (e.g. Malkus 1954 a), and this provides 
strong support for the concept of surface convection layers within which 
most of the temperature variation occurs. 

The concept of a surface layer with a structure determined by the heat 
flux through it is very similar to the concept of a constant stress layer in 
turbulent flow near a wall, and this suggests that there may be a part of the 
surface layer within which the direct effects of conductivity and viscosity 
are negligible. If there is such a region, within it we have 


as 


log = (3.7) 
Fa = Bo (3.8) 


where the constants may depend on Prandtl number. ‘The prediction 
(3.7) is supported by measurements of heat transfer in the lower atmosphere 
(Priestley 1954). 

The only complete theory of the convection between parallel planes is 
that proposed by Malkus (1954 b) and suggested to him by his observations 
of successive transitions between the initial cellular convective flow and the 
final state of completely disordered turbulent motion. Malkus suggests 
that the flow adjusts itself to transfer the maximum amount of heat 
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compatible with the boundary conditions and subject to the additional 
restrictions, 

(a) that the heat flux is everywhere down the gradient of mean 
temperature, 

(6) that the variations of temperature and vertical velocity may be 
represented by Fourier series which terminate after a finite number of 
terms. (‘The number of terms in the two Fourier series is determined by 
requiring that the distribution of mean temperature, which ts defined by 
these conditions, should be neutrally stable with respect to perturbations 
of the smallest wavelength occurring in the Fourier series.) 

Unlike most current theories of turbulent motion, this theory is based 
on generalizations more like to those of thermodynamics than those of 
fluid mechanics, but it has the considerable virtue of making definite 
numerical predictions. For our purpose the most relevant are (i) that, 
not too close to the bounding surfaces, 


1 dT 1 « flog 
or, more usefully, 
K T,\2 


where 77 = 7,7,, and (ii) that the mean square of the temperature 
fluctuations, averaged over the whole space, is 


2 T\3 « OD 


where the double bar signifies a double average, first over time at a point and 
then over ali values of z in the field. It should be remarked that the Malkus 
theory is consistent with the assumption of independent surface layers at 
high Rayleigh numbers but not with the additional assumption that the 
direct effects of viscosity and conductivity are negligible. By using 
published measurements of the heat transfer coefficient, which may be 
represented by equation (6.2), equation (3.10) may be written as 


T (2000v/x)!2 
logge = (2) (3.12) 


a 27” 


valid for values of z/sy within the surface layer which are not too small. 
This is of the same functional form as equation (3.4). 


4. EXPERIMENTAL ARRANGEMENTS 
The source of heat in these experiments was an electrically heated 
rectangular duraluminium plate of dimensions *30 cm x 40 cm x 1 cm. 
In order to secure a uniform distribution of surface temperature and to 
allow direct measurements of heat transfer, this plate formed the upper side 
of a duraluminium-asbestos-duraluminium sandwich (figure 1), the lower 
plate being heated by a grid of resistance wire wound on an asbestos former 
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and placed against its lower side. Four copper-constantan thermocouples 
were inserted in each duraluminium plate, and from the mean temperature 
ditference and the measured thermal conductivity of the asbestos thermal 
resistance the total heat loss from the top plate could be calculated. From 
this was subtracted the heat loss by radiation and by conduction along the 
pillars supporting the upper surface. The radiation correction was 
computed using the emissivity 0-24 appropriate to polished aluminium, 
and was found to form about 30°, of the total heat loss. In the parallel plate 
experiments, the upper boundary was a similar aluminium plate forming 
part of a water jacket. The temperature of this plate could also be measured 
by thermocouples set in it. Both plates were ground flat after assembly 
and were kept at known distances apart by glass pillars ground to size. 
Vertical walls of rubber or asbestos were always fitted to prevent entry of 
air from the edges of the space. For the measurements over a single heated 
surface, vertical side walls, 60 cm high, were provided for the same purpose. 
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Figure 1. Experimental arrangement : diagrammatic section. 


‘Temperature measurements in the air space were made using resistance 
thermometers of Wollaston wire, 0-00025 cm in diameter of measured 
temperature coefficient 0-00353 deg !, along which were passed measuring 
currents of about one milliampere. The length of the sensitive, etched 
portion of the wires was usually about 3 mm, giving a resistance of about 
30 ohms. ‘The response time of such a wire is of order one millisecond. 
The wires were soldered between two 30 s.w.g. copper wires stretched 
horizontally between supports on a movable carriage spanning the whole 
apparatus. ‘This carriage could be moved either vertically or horizontally 
by micrometer screws. 

‘The resistance of the Wollaston wires was measured with a Wheatstone 
bridge supplied with alternating current of 1000 ¢.p.s. The out-of-balance 
signal of the bridge was amplified and supplied to a phase-sensitive detector. 
‘The current output from such a detector is a linear function of the wire 
resistance and, passed through a reflecting galvanometer of period about 
one second, gave deflections of a light beam that could be recorded by a drum 
camera (see figure 2, plate 1 and iigure 9, plate 2 for some sample records), 
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The whole system was calibrated by substituting known resistances for the 
thermometer element. From these records, the mean temperature was 
determined using a planimeter to measure the area under the trace, and 
the variance by measuring the deviation at intervals along the record. 
The duration of each record was about eight minutes. 
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Figure 3. Distribution of mean temperature between parallel plates 
(D = 2:92'cem, = 52:4°¢, = 12°37 €, OSS 10). 


In general, the supports of the wire were not at the same height as the 
wire itself and so were not at the same temperature. Anappreciable amount 
of heat may be conducted along the supports and lead to a systematic error 
in mean temperature. A correction for this was calculated and applied to 
the measured values. 

Some measurements of the autocorrelation function of the temperature 
fluctuations were also made. For this, the rectified and smoothed output 
from the Wheatstone bridge was used to modulate the frequency of an 
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oscillator, and the oscillations recorded on magnetic tape. ‘This record 
could be played back and the autocorrelation function determined by 
comparing frequencies at points on the tape separated by a known fixed 
distance, that is, the frequencies recorded at moments separated in time by 
the time necessary for the tape to travel the fixed distance. The actual 
comparison consists of the observation of the rate of coincidences between 
standard duration pulses constructed from each set of played-back oscillations, 
the rate being proportional to the product of the instantaneous pulse 
frequencies (‘Thomas 1956). 
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Figure +. Distribution of mean temperature between parallel plates 
(D: ='4°88.cm, 7, 46-8° C, 7, 10°6" C, A X10"). 


5. ‘THE CONVECTION BETWEEN PARALLEL HORIZONTAL PLANES 


Three sets of measurements were made with very nearly the same 
temperature ditference between the plates and three different spacings. 
‘Table 1 summarizes the experimental conditions and the measured heat 
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Figure 9. ‘Temperature fluctuation records over a heated horizontal plate 


(T, = 316°C, Tp = 21°6°C). 


4 


Turbulent convection over a heated horizontal surface 481 


fluxes. It will be noticed that the measured coefficients of heat transfer 
are rather less than those found by Malkus (1954 a) who used water and 
acetone, but that the difference seems to decrease with increasing Rayleigh 
number. 

The measurements of mean temperature are shown in figures 3 to 5. 
It will be noticed that the profile at the highest Rayleigh number is markedly 
asymmetrical and even contains a section with reversed temperature gradient. 
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Figure 5. Distribution of mean temperature between parallel plates 
5:37 eu, = 465° €, = C, A — 6S x10). 


It is believed that this is due to the onset of a circulatory motion of dimensions 
similar to those of the whole apparatus, with air rising from the centre of the 
lower plate where the measurements were taken. This motion would cause 
a finite mean shearing stress at the boundaries and logarithmic distributions 
of mean velocity and temperature may be expected. Ifthe mean temperature 
is plotted against the logarithms of the distances from the upper and lower 
surfaces (figure 6), substantial regions of logarithmic variation are found. 
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If the other two profiles are plotted in like manner, there is no sign of a linear 
section on the curves. ‘The onset of this circulation at a Rayleigh number 
around 5-10° is probably due to the side walls forcing the large eddies of the 
convection to remain relatively fixed in position. It is probable that, 
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Figure 6. Logarithmic temperature profiles (A = 6°75 = 10°). 


for greater ratios of plate width to separation, thesé eddies could wander 
freely and then measurements at a fixed point would be identical with 
spatial means—which clearly they are not, for this particular profile. 
For comparison with the Malkus theory, the three profiles have been 
plotted against cot(7z/D) and straight lines fitted to the central portion 
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(figure 3 to 5). At the lower Rayleigh number there is a substantial linear 
region but the slope is less than the predicted value. ‘The slope of the 
second profile agrees more closely with the prediction, but substantial 
departures from linearity exist. As shown by this method of plotting, these 
departures are similar in kind to the gross distortion in the third profile 
and may be caused by a weak circulation of large scale. 
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Figure 7. Temperature fluctuation intensities between parallel plates 
(4 = 6°75 x 105). 


Measurements of the root-mean-squares of the temperature fluctuations 
for the highest Rayleigh number are shown in figure 7. The average level 
is 1:9 C, compared with the value of 49° C given by equation (3.11). 
This discrepancy is probably too large to be explained by the finite time of 
response of the recording equipment. As would be expected from the shape 
of the mean temperature profile, the intensity of the temperature fluctuations 
reaches a peak just outside the conduction layer. 
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Figure 8 shows the autocorrelation function of the temperature fluctua- 
tions, defined as 
R(r) = O(t)A(t + (5.1) 
at three values of z. It will be noticed that the time scale of the fluctuations 
decreases sharply as the lower surface is approached. 
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Figure 8. Autocorrelation functions of the temperature fluctuations 


(D = 5-87 cm, T,—T, = 38-8°C). 


6. "THE CONVECTION OVER A SINGLE HORIZONTAL PLANE 

Before discussing these measurements, it is necessary to consider the 
relation of the experimental arrangement—a horizontal heated plane forming 
the bottom of an open box of roughly cubical form—to the infinite heated 
plane considered in the theory. Unless there is appreciable large-scale 
circulation, and subsidiary experiments showed that there were not, the 
conditions near the surface will approximate to those near an infinite plane 
if the distance from the surface is small compared with distance from the 
walls and distance from the top of the box. With the single exception of the 
autocorrelation measurements, no measurements were made more than three 
centimetres from the surface and the condition is satisfied. 


Measured quantities Derived quantities 


Hx10* Hx10- T, | log7,/T, Q1076,x10? 2, 


316°C 21:6°C 0°535 -21-6° C} 0-0434 1-48 1:12 0-158 
236°C 1:79 24-4°C| 0:0917 499 2-80 0-116 
3:38 3:09 28-5°C| 0-1565 9:20 4:44 0-100 


Table 3. Convection over a horizontal plane’ (c.g.s. units). 


The boundary conditions of the theoretical flow are defined by the 
temperature of the plane and the temperature at infinity, the ‘ambient 
temperature’, which corresponds with the temperature at the base of the 
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convection plume that rises from the top of the box and carries away the heat 
convected from the plane. Unfortunately, measurements of this base 
temperature, which may differ from room temperature by as much as 
8’ C, were not made at the time of the experiments, and it is necessary to 
estimate it by extrapolation of the measured temperature profiles. In 
figure 10, the three measured temperature profiles are plotted non- 
dimensionally as #)'log(7/)T,) vs 3/39, using the room temperature as 
reference temperature. It is clear that vertical displacement of the curves, 
corresponding to changes in the effective ambient temperature, will bring 
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Figure 10. Distribution of mean temperature over a heated horizontal plate 
(similarity plot referred to room temperature). 


them into near coincidence. ‘lhe results plotted in figure 11 have been 
adjusted in this way, using the smallest possible positive differences of 
ambient temperature from room temperature, that is, zero for the first 
profile for which 7,— 7, = 10-0 C. The necessary adjustment increases 
rapidly with plate temperature, probably because air currents in the room 
are more effective in ventilating the top of the box when the convection 
plume is weak. 

The universal non-dimensional profile so obtained may be represented 
within experimental error by 


651 log(T/T,,) = (6.1) 


for all values of s/s, greater than two (see figures 11 and 12). It is not 
possible to represent the profile by a s-'3-law over a substantial range of 
:/3, Without assuming ambient temperatures much less than room 
temperature and an upper limit to its validity at s/s, = 12. ‘This upper 
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limit w 
dominant at only 1 cm from the plate, which seems highly unlikely, and 
that the effective ambient temperature is less than room temperature, 
which 1s equally unlikely. 
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ould imply that the influence of the distant boundaries becomes 
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plot referred to a selected ambient temperature). 
equation (6.1). 
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Figure 12. Mean temperatures plotted against (z/z)~! and (z/z)-!/° 
for z/z9 > 2. 
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Direct measurements of heat transfer are shown by plotting first 

#5 log(7,/Tp) and 67! log(7,/T,,) against log(7,/7T,) (figure 13). In the 

second graph, appropriate values of 7,— 7, have been obtained by inter- 

polating between the values used in figure 11. If the convection is really 

independent of the side walls and the distant boundaries, 67! log(7)/7,,) 

should be independent of temperature difference, as is implied by equation 

(3.3). With these values for the ambient temperature, these measurements 
indicate that 

3\ 14 
log = 3-49, = 3-4 (=) (6.2) 
a KS 
nearly independent of temperature difference. ‘The coefhicient 3-4 may 
be compared with the value 3-1 obtained from the temperature profiles. 
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Figure 13. Measurements of heat transfer from a heated horizontal plate (plotted 
non-dimensionally). 


It is also interesting to compare these values of heat transfer with measure- 
ments between parallel planes at large Rayleigh numbers. ‘The 
measurements of Malkus (1954 a) and others are closely represented by 
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If the heat transfer at each surface were little affected by the presence of 
the other surface, we would expect the ratio 
log (T,/T,,) 
(Q°, kg)! 4 
to be a constant and the temperature midway between the surfaces to be 
the geometric mean of the surface temperatures. From equation (6.3) 
and inserting the value of the Prandtl number v/« for air, we find 
T 3\ 1/4 
log =" = 3-12 2 (6.4) 
= (7, in good agreement with the previous values*. Even 
the measurements of heat transfer of the previous section give values of 
4-0, 3-7 and 3-6 at Rayleigh numbers of 0-85 x 10°, 3-77 « 10° and 6-75 x 10°, 
although these are too low for independence of the two surface layers. 
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Figure 14. Temperature fluctuation intensities over a heated horizontal plate. 


Measurements of root-mean-square temperature fluctuation, plotted 
non-dimensionally as against (gQ/«*)!4z, are shown in 
figure 13. ‘The measurements for the two heat transfers only agree at the 
larger values of (gO/«*)!4z. It is very probable that the response time of 
the galvanometer prevented the complete recording of the rapid fluctuations 
occurring near the surface. The single measurement of the autocorrelation 
function (figure 14) shows that a substantial part of the intensity resides in 


* These non-dimensional equations refer to air with v/k = 0°77. Accepting 
equation (6.3) as correct, the general form of equation (6.4) is 
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Fourier components of periods less than the period of the galvanometer, 
one second. It may be mentioned that more recent work using equipment 
with faster response has registered temperature fluctuations of much 
greater amplitude. 
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Figure 15. Autocorrelation function of the temperature fluctuations over a horizontal 
plate (7,—Tp = 70-9°C, z = 4 cm). 


7. DISCUSSION OF RESULTS 

‘The most interesting feature of these measurements is the close 
correlation between heat transfer from a single surface and heat transfer 
between parallel planes. Even though the Rayleigh numbers of the parallel 
plane experiments were comparatively small, they tend to confirm the 
experimental finding of Malkus that, at high Rayleigh numbers, the heat 
transfer coefficient is nearly proportional to the cube root of the Rayleigh 
number and so is independent of the separation of the surfaces. This 
must mean that the conditions near each surface dictate the heat transfer 
and that the central region is a region of negligible temperature gradient, 
easily capable of transmitting the imposed heat flux. Nearly all the 
temperature drop occurs in the wall layers which have a temperature 
distribution determined by the heat transfer, the conductivity (and viscosity) 
and the gravitational field. Additional confirmation of this view is found 
in the observation of Malkus (1954 a) that the heat transfer between parallel 
planes is practically unaffected by the presence of vertical partitions even 
though they restrict greatly the motion of the central region. 

Although both the existing theories of convection conform to this 
general pattern, neither is capable of describing accurately the measured 
profiles. [he similarity theory (Priestley 1954) requires that log(7/7T,) 
should be proportional to s-'* for al/ large values of 3/2, where 7, the 
temperature at large distances from the surface, is certainly not less than room 
temperature. It is impossible to fit any substantial part of the observations 
by such a power law without using a value of 7), less than room temperature. 
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‘The very interesting theory of Malkus (1954 b) is more successful in as 
much as its prediction that 


log = (7.1) 
A 2 
fits the observations well if C = 3-0. However, consistency of the heat 
transfer measurements and the theory requires that C = 2-0 (see equation 
(3.9)). 

The failure of the similarity theory to describe the temperature variation 
within the range of measurement implies that the local temperature gradients 
and other local characteristics of the convection are affected either by the 
conductivity (or viscosity) of the fluid or by the convection near the distant 
boundary. ‘This may appear surprising in view of the considerable success 
of the similarity treatment of flow in the constant stress region of turbulent 
wall flow, for which similar assumptions of independence of viscosity and 
presence of distant boundaries are made. ‘The analogy fails in that the 
constant stress layer is a layer of substantially uniform turbulent intensity 
in which no appreciable net transport of turbulent energy is expected, while 
the surface layer, as described by the similarity theory, would have a turbulent 
intensity varying as (Ogz)?* and substantial transport of turbulent energy 
towards the surface would be expected. Also, the intensity of the tem- 
perature fluctuations varies considerably and they will be transported through 
the surface layer. ‘These considerations (and the present results) suggest 
that interaction between ditferent parts of the surface layer may be more 
effective than the similarity theory supposes. It is possible that the 
conditions assumed may exist for values of z/s) beyond the range of measure- 
ment, where the difference from ambient temperature is very small compared 
with the temperature difference across the whole layer. If the 
meteorological evidence is accepted, it presumably means that these 
measurements were made in such a region. 

The defects in the Malkus theory are difficult to diagnose, partly because 
the theory is based on general statements not easily related to the more 
common approach through the equations of motion and heat. ‘The difference 
between the theoretical and observed distributions of temperature in the 
surface layer indicates that the convection has not attained the optimum 
condition specified in the theory. In fact, it seems to have settled down 
with a heat transport about 70°,, of the maximum set by the theory. ‘Two 
possible causes are (a) that the requirement that the minimum temperature 
gradient is zero is insufficiently restrictive, and (4) that the Fourier series 
used to describe the convection does not satisfy all the boundary conditions. 
In the central region at the higher Rayleigh numbers, there is a systematic 
tendency to lower gradients, which is a departure in the opposite direction. 
Some of this may be due to the presence of a large-scale circulation, but 
measurements of the closely analogous flow between concentric rotating 
cylinders show an extensive region of constant angular momentum, the 
quantity analogous to temperature in convection (‘Taylor 1935). Pai (1943) 
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finds that the gradient of angular momentum is reversed in two small regions 
adjacent to each cylinder, and explains this as a result of large-scale circula- 
tions which are part of the general flow. It may well be that the large-scale 
circulation postulated to explain the form of the temperature profile at the 
largest Rayleigh number is typical of the flow as a whole, except in its 
position, which is determined by the presence of side walls. For the 
present, it may be concluded that the Malkus theory is in fair but not good 
agreement with experiment, which is a remarkable achievement for a theory 
of turbulence not containing any disposable constants. 


‘The experimental work described above was dene in the Cavendish 
Laboratory by the first author (D. B. ‘T’.), who is pleased to acknowledge 
the sponsorship of the work by the Meteorological Office. ‘lhe second 
author has contributed some additional analysis of the results. 
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SUMMARY 

In an earlier paper the uniform shear flow past a sphere was 
studied, by investigating how vortex lines are deformed by the 
‘primary flow’ (flow in the absence of shear), and deducing the 
‘secondary’ vorticity field (first approximation for small shear). 
In another paper the image system associated with each element 
of secondary vorticity was found, whence the Biot-Savart law 
can be used to determine the secondary flow field by integration. 
The integration is here carried out for the ‘downwash’ (secondary 
flow component perpendicular to the undisturbed flow, down the 
velocity gradient) on the dividing streamline. Difficulties due to 
the infinite domain of integration and singularities of the integrand 
are overcome by selecting variables of integration carefully and 
using known analytical properties of the secondary vorticity. 
From the computation of downwash is inferred the first approxi- 
mation (for small shear 4) to the ‘ displacement’ 6 (displacement 
of the dividing streamline, up the velocity gradient, far upstream 
of the sphere). If U is the upstream flow velocity and a the radius 
of the sphere, the computed value of lim(U/6/Aa?) is 0-9. 

A—0 

Details of the calculation show that the secondary trailing 
vorticity is not an important contributor to the displacement. 
‘The downwash is due almost entirely to vorticity upstream of the 
sphere (Hall’s earlier simplified theory gave good results, e.g., 1-24 
instead of 0-9, because it concentrated on the effect of local vorticity 
in producing downwash); further, this produces displacement 
principally through its image vorticity. 

The relation between theories for a sphere and experimental 
results on Pitot tubes (beginning with Young & Maas 1936) is 
discussed. Theoretical evidence on tertiary- and quartary-flow 
effects is used here in the light of recent work which renders the 
successive-approximation sequence uniformly valid at infinity. 
The conclusion is that the theories, taken together, are not incon- 
sistent with the experimental evidence that (i) at values of the 
‘shear parameter’ da U at which the displacement is measurable 
the ratio 6/a seems to have asymptoted to an approximately constant 
value, and (ii) displacement is greatly reduced in supersonic flow 
(Johannesen & Mair 1952) or when ‘sharp-lipped’ tubes are 
used (Livesey 1956), 
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finds that the gradient of angular momentum is reversed 1n two small regions 
adjacent to each cylinder, and explains this as a result of large-scale circula- 
tions which are part of the general flow. It may well be that the large-scale 
circulation postulated to explain the form of the temperature profile at the 
largest Rayleigh number is typical of the flow as a whole, except in its 
position, which is determined by the presence of side walls. For the 
present, it may be concluded that the Malkus theory is in fair but not good 
agreement with experiment, which is a remarkable achievement for a theory 
of turbulence not containing any disposable constants. 


‘The experimental work described above was done in the Cavendish 
Laboratory by the first author (D. B. ‘T’.), who is pleased to acknowledge 
the sponsorship of the work by the Meteorological Office. ‘lhe second 
author has contributed some additional analysis of the results. 
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SUMMARY 

In an earlier paper the uniform shear flow past a sphere was 
studied, by investigating how vortex lines are deformed by the 
‘primary flow’ (flow in the absence of shear), and deducing the 
‘secondary’ vorticity field (first approximation for small shear). 
In another paper the image system associated with each element 
of secondary vorticity was found, whence the Biot-Savart law 
can be used to determine the secondary flow field by integration. 
The integration is here carried out for the ‘downwash’ (secondary 
flow component perpendicular to the undisturbed flow, down the 
velocity gradient) on the dividing streamline. Difficulties due to 
the infinite domain of integration and singularities of the integrand 
are overcome by selecting variables of integration carefully and 
using known analytical properties of the secondary vorticity. 
From the computation of downwash is inferred the first approxi- 
mation (for small shear A) to the ‘ displacement’ 6 (displacement 
of the dividing streamline, up the velocity gradient, far upstream 
of the sphere). If U is the upstream flow velocity and a the radius 
of the sphere, the computed value of lim(U/6/Aa?) is 0-9. 

Details of the calculation show that the secondary trailing 
vorticity is not an important contributor to the displacement. 
‘The downwash is due almost entirely to vorticity upstream of the 
sphere (Hall’s earlier simplified theory gave good results, e.g., 1-24 
instead of 0-9, because it concentrated on the effect of local vorticity 
in producing downwash); further, this produces displacement 
principally through its image vorticity. 

The relation between theories for a sphere and experimental 
results on Pitot tubes (beginning with Young & Maas 1936) is 
discussed. ‘Theoretical evidence on tertiary- and quartary-flow 
effects is used here in the light of recent work which renders the 
successive-approximation sequence uniformly valid at infinity. 
The conclusion is that the theories, taken together, are not incon- 
sistent with the experimental evidence that (i) at values of the 
‘shear parameter’ Aa/U at which the displacement is measurable 
the ratio 6/a seems to have asymptoted to an approximately constant 
value, and (ii) displacement is greatly reduced in supersonic flow 
(Johannesen & Mair 1952) or when ‘sharp-lipped’ tubes are 
used (Livesey 1956). 
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1. INTRODUCTION 

In this paper the mechanism of the ‘ Pitot-tube displacement effect’ 
is analysed by further study of the shear flow past a sphere, following on 
three earlier papers, ‘* The image system of a vortex element in a rigid sphere” 
(Lighthill 1956 a, to be referred to hereafter as I), “* Drift’? (Lighthill 1956 b, 
with corrigenda in Lighthill 1957 a, to be referred to collectively as D), 
and ‘The displacement effect of a sphere in a two-dimensional shear flow” 
(Hall 1956, to be referred to as H). 

When a parallel shear flow, such as a boundary layer or wake, is 
investigated with a Pitot tube (that is, a tube with the open end pointing 
upstream and the other end closed by a manometer), it is found that the 
pressure in the tube is greater than the total pressure on the streamline 
approaching the centre of the tube orifice; it is equal rather to the total 
pressure on a streamline displaced in the direction of higher velocities, by 
an amount 6 usually known as the ‘displacement of the effective centre’ 
of the Pitot tube. For the tubes whose ratio of internal diameter d; to 
external diameter d, is 0-6, Young & Maas (1936) found that 6/d, was 
scattered fairly randomly about 0-18, and Macmillan (1956) obtained an 
average value 0-15 for 6/d.. 

Two recently observed phenomena enhance the desirability of 
understanding the displacement effect physically. First, in supersonic 
flow, Johannesen & Mair (1952) showed that the displacement effect 
practically vanished in a wake for MW = 1-96. Understanding of this result 
is made vastly harder by the presence of the ‘ bow shock wave’ in front of the 
Pitot tube, but in view of the considerable region of subsonic flow between 
that shock wave and the orifice, it is safe to say that no explanation can 
reasonably be attempted until the incompressible-flow phenomenon is 
understood. 

Secondly, although Young & Maas (1936) found some indications that 
6/d, increased slightly with increase of d,/d,, Livesey (1956) has made 
experiments which indicate that there is practically no displacement when 
d;=d,. This condition was obtained by chamfering the front part of the 
tube wall to a cone-frustum shape, so that the orifice became sharp-edged. 
As a check on his experimental technique, Livesey obtained results for 
d,d, = (0-6 comparable with those noted above. The contrasted results 
are difficult to understand without a rather full discussion of the mechanism 
of the displacement effect. 

It might be thought that useful data for such a discussion would be 
found by a study of the analogous two-dimensional problem, in which a 
two-dimensional ‘ Pitot channel’, whose external thickness becomes 2c 
downstream of the nose, is placed in a parallel shear flow. However Hall 
shows (H, Appendix) that the displacement of the dividing streamline in 
this problem is given by 
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where U + Ay is the velocity distribution in the oncoming stream*. For 
values of U/Ac occurring in practice, this result is approximately 


Ad, 


Ww 

a value well below the experimentally determined values of the displacement 
for the values Ad,/U — } which are used in practice. Similar results are 
found for two-dimensional shear flow about cylinders of various shapes 
(see Mitchell & Murray 1955 as well as references in H). Hall shows, 
however, that a three-dimensional theory, which takes into account the 
stretching of the vortex lines of the oncoming shear flow as they pass over 
the body, but still neglects viscous effects, can be used to obtain displacements 
of the right order of magnitude. ‘There is a simplifying approximation in 
Hall’s work, which will be critically examined below by comparison with 
a theory which avoids the approximation, but the results of this examination 
will be found to be favourable. 

Before further detailed discussion it may be worth indicating where the 
subject stands within hydrodynamics as a whole. It falls under the general 
heading ‘secondary flow’ which is used to describe what happens when a 
parallel shear flow is disturbed in a three-dimensional manner (for example 
by being led round a bend in a pipe or channel; or deflected by a pressure 
gradient transverse to the streamlines; or subjected to Coriolis force as 
in meteorology ; or confronted with an obstacle as in the present problem). 
More particularly, it is used to denote the departure of the disturbed flow 
from a so-called ‘ primary flow’, in which the streamlines are the same as 
if the original parallel flow had been uniform instead of sheared. 

There are two useful ways of looking at secondary flows. In the older 
approach (see, for example, Goldstein 1938, chap. 2) one considers what 
unbalanced pressure gradients are implied by the primary flow when 
account is taken of inertial forces, altered in magnitude because the original 
parallel flow is not really uniform. ‘The secondary flow is regarded as 
maintained by such pressure gradients, and limited by viscous resistance. 
This view of the matter, though always of qualitative value, leads to useful 
quantitative results only at rather low Reynolds numbers, when the inertial 
forces may be regarded as a small perturbation of a régime dominated by 
viscous resistance (Dean 1927, 1928; Cuming 1951). 

In flows at high Reynolds numbers it has therefore become usual 
(Squire & Winter 1951; Hawthorne 1951, 1954; Hawthorne & Martin 
1955; Preston 1954) to seek quantitative results by another method, 
namely study of the vorticity field. ‘The technique used can be interpreted 


*In the last sentence of this Appendix, note that “‘on the body” should 
read as ‘‘at infinity’’. (To make % = O the dividing streamline, 4%, must tend to 
—4Ac? as the external body surface becomes flat. But a bounded harmonic function 
must tend to the same limit at infinity in all directions. Hence also upstream 
—4Ac*, so that ~ 3Ay?—5Ac*® and the displacement 6 of the dividing 
streamline ys = 0 satisfies (A.9).) Also, in equation (A.8), $4 should read as }Ay’. 
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geometrically as a study of how vortex elements in the original shear flow 
are stretched and rotated by the primary flow, yielding a ‘secondary vorticity 
tield’, whose associated velocity field is the ‘secondary flow’. In D, a 
mathematical approach was described which, while essentially equivalent 
to those of the authors cited, corresponds rather more directly to this 
geometrical picture. 

Such a method is obviously of only approximate validity, partly because 
it neglects diffusion of vorticity by viscosity or turbulence, and partly because 
the vortex elements are in reality stretched not by the hypothetical primary 
How but by the exact flow. The latter inaccuracy can in principle be reduced 
by calculating a ‘tertiary vorticity field’ by considering the stretching and 
rotation of vortex lines by the combined primary and secondary flows, 
and deducing a ‘tertiary flow’ from it, and one can even consider quartary 
Hows, as Hall was able to do (H, p. 154) by the use of his simplifying 
approximation. 

The vortex approach to secondary flow has so far had the greatest 
success with internal flows (through bends in pipes, or cascades), for which 
the calculation of the velocity field from the vorticity field can be fairly 
straightforward. With external secondary flows, such as those resulting 
from parallel shear flow past an obstacle, this calculation is more difficult, 
and work apart from Hall’s (H) has been limited to calculations of the 
secondary vorticity field (Hawthorne 1954; Hawthorne & Martin 1955; 
and D, §2 and §7). ‘The problem of how to deduce the secondary velocity 
field without making Hall’s approximation is considered in some detail in 
this paper, partly to establish the value of Hall's approximation, partly 
to improve understanding of the effect, and partly to improve techniques 
for the mathematical calculation of three-dimensional, fully rotational 
flows. 

The calculations are still confined to the case of a sphere, which can be 
regarded in many ways as a typical bluff obstacle. ‘The work was consciously 
prepared for in D by the calculation of the secondary vorticity field, and 
in I by the calculation of the image system of a single vortex element in the 
sphere. ‘This type of image approach, combined with the Biot—Savart law, 
is found to be more convenient in th. present problem than any approach 
based on images of complete vortex lines. 

The relevance of calculations for a sphere to the Pitot-tube problem is 
supported by the finding that vorticity downstream of the sphere contributes 
negligibly to the displacement effect. ‘This also kills a plausible hypothesis, 
namely that secondary trailing vorticity might be largely responsible for the 
displacement effect, which therefore would be greatly reduced at supersonic 
speeds for which the trailing vorticity has no upstream influence. 

A more correct physical interpretation of the results, as discussed in § 4, 
is that in all cases the image system of the vorticity upstream of the obstacle 
is the main cause of the displacement. It is shown in $5 how this can be 
regarded as responsible for the altered situation when either a supersonic 
main stream or a sharp-lipped orifice is used, 
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2. EVALUATION OF THE SECONDARY DOWNWASH AHEAD OF THE SPHERE 

‘The shear flow past a sphere will here be studied in the notation set out 
in D. ‘hus, Cartesian axes are used, with origin at the centre of the sphere, 
such that far upstream the velocity field is 

v, = U+Ay, v, = v, = 0, (1) 
where suffixes denote components. In addition, spherical polar coordinates 
r, #, A such that 

x =rcosé, y = rsin@cosa, z= rsin@sinaA (2) 
are used; in these coordinates the sphere is r = a. 

The primary flow is the irrotational flow about the sphere, associated 
with an upstream velocity field as in (1) but with 4 = 0. ‘The secondary 
flow represents the first-order correction for non-zero A: thus, it is the 
term in Aa U ii an expansion in powers of that non-dimensional parameter. 

The secondary vorticity field was determined in D, §7. Certain 
asymptotic properties of the secondary velocity field were derived in D, § 3. 
Here, such further properties of it are derived as are necessary to determine 
the displacement of the stagnation streamline to the first order in the 
parameter Aa/U. 

For this one must know the distribution along the line y = z= 0, 
x < —a (which is the stagnation streamline for A = 0) of the secondary 
flow component normal to that line. ‘This component is actually in the 
negative y-direction. It may be referred to as ‘downwash’ if we think of 
the upstream velocity distribution as horizontal and increasing upwards. 
Accordingly, we write D(s) for the value of (—7,) at the point (—s, 0, 0), 
and call D(s) the ‘ downwash function’. 

Now, the secondary flow may be calculated in three parts as explained 
in D, $3. First, there is a part, which is the gradient of a potential ¢, 
vanishing at infinity, and whose normal velocity component on the surface 
just cancels out that associated with the uniform shearing motion 7, = Ay, 
v, =v, =. In the case of a sphere this condition may be satisfied by a 
single solid harmonic (easily obtained by inspection), namely 

Aa’xy 
(3) 
This part may be regarded as the velocity field due to the image system of 
the undisturbed vorticity field (0,0, —A) in the sphere r=a; we see 
that this image system consists of a single quadrupole at the centre. ‘The 
contribution of this part to the downwash function D(s) is 
Ds) = (4) 

Secondly, there is the Biot-Savart field of the,‘ vorticity change’ o,. 
that is, the difference between the secondary vorticity field and its undisturbed 
value (0,0, —A). The contribution of this Biot-Savart field (D, (82)) to 
D(s) is 

1 W1,%— +5) 
| +y? + 27}? 
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dxdydz. (5) 
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‘Yo evaluate (5) it is convenient to change over to spherical polar coordinates, 
when it becomes 


w,,8sin dA + ay cos r)sin A + +rcos 8)cos 
(s* + 2sr cos + r?)3* 


x (6) 


‘Thirdly, there is a further irrotational part of the secondary velocity 
field, which must be added to the Biot-Savart field of w, so that together 
they may satisfy the boundary condition on the surface of the sphere. 
In I this was shown to be simply the Biot-Savart field of a certain system 
of image vorticity inside the sphere. Using the term ‘strength’ to denote 
the product of the volume dl’ of an elementary region with the vorticity 
inside it, it was shown that transverse vortex elements of strengths w,dV 
and w,dl’ have images of strengths (—a/r)w,dV and (—a/r)w,dV  res- 
pectively at the inverse point (a?/r,0,A) while a radial vortex element of 
strength w,dV has an image system consisting of (i) a radial vortex element 
of strength (+a r)w,dV at the inverse point, and (ii) a uniform line vortex, 
of strength (—1l/a)w,dV per unit length, stretching between the inverse 
point and the origin. 

It follows by comparison with (6) that the contribution of this image 
vorticity to the downwash function D(s) 1s 


D.{s) | | | 


cos + a? 'r)sinA— w,;\s + (a? r)cos A 


s* + 2s(a?/r)cos + (a?/r)?)3? 


sina 


as sin + 2s(a? r)cos + (a? r)*} 
(7) 


where the coefficient of (—w,,/a) inside the square brackets is obtained by 
integrating with respect to r from 0 to a* r (the value of 7 at the inverse 
point) the coefficient of @,, 1n the curly brackets in (6). 

Now, @,, takes the simple form (D, (1)) 


( 3\ -12) 
= (A cos), 1 (1 = (8) 


Asa result, the terms involving w,; in D,(s) + D,(s), which may be designated 
as 1),(s), can be greatly simplified. Since 
(s +rcos#)sin 6 dé (s >r) 
) 


(9) 


we have 


1 3\-12) q°\ -12) 
D,(s) = al | (1 - dr — | <|- (1 3) rar ar | 


s\32) 
a 


J 
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the last term (independent of s) being obtained by evaluating the infinite 
integral in terms of factorial functions. 
The vorticity components w,, and w,, can be written 


w,, = (Asind)F,(r, 4), wy = (AsindA)F,(r, (11) 


where in D, $7 the quantities F, and F, were tabulated for points on three 
particular streamlines, and determined asymptotically near the axis of 
symmetry, near the surface of the sphere and far from the sphere. In 
terms of F, and F,, the terms involving «,, and wy, in D,(s)+ D,(s), which 
may be designated as D,(s) and D,(s) respectively, may be written as follows: 


D,(s) = 44 


~a 


ssin@ 


(s? + 2sr cos 6 + r?)8? 


a ssiné | 
+ 2s(a*/r)cos 6+ assin®é 


scos 6+ a?/r 
x Hie — 12 
( ) 


scosO+r 


a scos0+a?/r 
+ 2s(a*/r)cos 6 + at/r?}3* 


In the numerical evaluation of the integrals (12) and (13) difficulties 
arise from two main sources. First, there are the ordinary difficulties 
involved in the numerical evaluation of any integral when the domain of 
integration is infinite. ‘These (as we shall see) can be overcome in the 
ordinary manner, because of our rather full knowledge of the behaviour of 
the integrands for large r/a. 

The second source of difficulties is the unboundedness of the integrands. 
For every s there is one value of r (namely r = s) for which one of the 
denominators in (12) and (13) vanishes on 6 = 7. In addition, F’, and F, 
have singularities, both on r = a, where F, is actually infinite, and on 6 = 0, 
where F, is actually infinite. Thus, there are difficulties in the neighbourhood 
of the whole dividing streamline (@ = 7, r = a, 6 = 0). 

This indicates the advantage of a change in the variables of integration. 
Instead of r, 6 we use py, 4, where po, introduced in D, §5 and defined by 
the equation (D, (60)) 


=) sin?4 = p>, (14) 


is constant along each streamline*. ‘hen the integration with respect to @ 
can be carried out first (and here the fact that thi$ integration is over a 
finite range () —< @ < w is an advantage); the singular character of the 
integrands comes in only when the second integration (with respect to pg) 


* Stokes’s stream function }Up% would do almost as well as ps as a variable of 
integration, but there are some small advantages in favour of py in the present problem. 
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is carried out, and there it is confined to the single point py = 0 which 
represents the whole of the dividing streamline. 

‘lo change the integrals (12) and (13) into integrals with respect to pg 
and 4, one has only to replace the differential product drdé in each by 


er r(r?— a? 
drd# = (— dp, dé = (15) 
‘The ranges of integration become 0 < 6 < 7,0 < py < ». 


The integrations with respect to @, keeping py constant, were carried 
out by Simpson’s rule, with equal intervals of 10° for 6, for three values 
(0-25, 0-5 and 1-0) of p)/a, using the values of r/a, F, and Fy, given in D, 
table 3, for the required combinations of values of @ and p,/a. At a few 
places, where the values of the integrands made it appear probable that 
Simpson’s rule was being slightly strained, intermediate values were 
computed and used to improve the approximation. 

If the result of such integration is written as D,(s, py) in the case of the 
integrand of D,(s), so that 


= | DAs, py) Apo, (16) 
~ 
and similarly with D,, the problem remains of carrying out the integration 
with respect to py. ‘Vhis was done by using Simpson’s rule in different 
forms for two parts of the integral, writing 


DAs, py) D,(s, py) dpy + | Ds, Py) Pa d(p; ') 


J 


l 
al D{s,0)+4D,(s, ja) + D,(s, + 


Ds 'a)+4a?D(s,a)+ lim'p; (17) 
3a : 
In this method one requires only those three values of D,(s, py) which 
were computed as described above, and the two limiting values D,(s, 0) and 
lim{p? D,(s,p))}.. An equation similar to (17) is used also for Dg. 


Now it may be shown that both D,(s,0) and D,(s,0) are zero. This 


result is deduced from the detailed behaviour of F, and Fy, as py > 0 
derived in D, $7 (equations (78) to (80)). The singularity at @= 7, r =s 
due to inverse (3 2)th powers in (12) and (13) produces no effect in D,(s, 0) 
because F’, vanishes on 6 = z, and none in D,(s,0) because of the factor 
scos@+r in the numerator. The singularity of F, on r = a 1s cancelled 
out by the vanishing on r = a of the term in square brackets in (13). The 
contribution to D,(s,0) as well as to D,(s,0) of this part of the streamline 
py = 0 therefore vanishes, since the coefficient of dp, dé in the differential 
element (15) vanishes as p, — 0 for fixed 4. Finally, the singularity in F, 
as @-+0 produces no effect in D,(s,0) because of the vanishing as 6 > 0) 
of the term in square brackets in (12). 


| 


Theory of the Pitot-tube displacement effect 501 
‘The limits as py ©,o0n the other hand, donotboth vanish. By(D, (77)) 


F, + ~ pycosecd, (18) 


r 


as py» &. Hence by (13) and (15) 
3 
Dils, po) ~ { - 3(<) sin?4 cos cosec*4 au cos 
“\Po 


but 
D,(s, py) ~ 4A | sin Od: O (20) 
so that 7Aa! 
lim{ p> po)} = Cae? lim {p> py)} = 0. (21) 


Por 


Note that the distant transverse vorticity is effective entirely through its 
images in the sphere (that is, only the second term in square brackets in 
(13) contributes to the asymptotic result (19)); the distant radial vorticity 
is less effective* because its image system consists of a nearly cancelling 
pair of vortex elements of opposite sense. 

Everything has now been found to enable D,(s) and D,(s) to be determined 
from equation (17). ‘This has been done for s/a = 1 and v2, and the results 
are given in table 1, together with values of D, (see (10)) and D, (see (4)), 
the contributions from the ‘ring vorticity’ w,, and from the image system 
of the undisturbed vorticity distribution. ‘lhe four terms are added to 
produce the total downwash given in the last column. 


| | | 
| Aa | D, Aa | D,/Aa | D,/: ‘Total D/da | 
} 

0-12, | 0-08, | 0-425 | 0-333 | 0-97 

y2| O11, | 0-04, | 0-083 | 0-083 | 0-33 | 

| 

| \ | 


Table 1 


For comparison with these values of D(s) we have for large s/a, by 


(D, (80), 
5 fa\2 
D(s) (—,,), sy =2=0 Aa, (22) 


where we have used the facts that the volume V,, of the sphere ts 47a’, 
and the volume V’, of fluid, whose mass is the ‘hydrodynamic mass’ 
associated with the sphere’s motion, is }V,. ‘Table 1 shows that (s/a)?D(s) 
is 0-97 for s/a = 1 and 0-66 for s/a = v2, so that it is tending fairly rapidly 
to the asymptotic value 2 = 0-42 as s/a—> x. 


* Actually, a more detailed analysis than that given above shows that 
DAs, po) = 


33 
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3. FIRST-ORDER DISPLACEMENT OF THE STAGNATION STREAMLINE 
‘The stagnation streamline, or ‘ dividing streamline’, must be in the plane 
of symmetry = = 0 (or A = 0) and on it 


v.rdd—vydr = 0. (23) 


For the irrotational flow about a sphere this equation becomes 


3 3 
cos dé + U sin + dr = 0, (24) 


which can be put into the form 


3 2 
rsin of \ 0, (25) 


4 


after multiplying through by U-'(1—a*/r*)'*. ‘The solution of (25) that 
must be used is 
rsiné{ 3 = (), (26) 


since other solutions give rsin@ —- « as r >a, which is impossible if the 
streamline is to meet the body. Hence # = 0 or z off the sphere r = a. 

When 4 40 we can obtain the displacement of the upstream part 
#=7 of the stagnation streamline by modifying the v7, term in (23) by 
including the downwash function D(s) of $2 init. The term in z, need not 
have a secondary flow term included, however, since it is multiplied by 
rdé which is itself a small* quantity of order 4 (vanishing, as we have seen, 
for 4d = 0). 

‘Thus, (24) is replaced by 


3 


Equation (27) after multiplication by U~1(1 — a*/r*)-1?, an integrating factor 
already used above, becomes 


a\ 12 D(r) dr 
( = 28 
OF (1 =) } (28) 


\ 
Hence the dividing streamline takes the form 


r D(s) ds 
12 [ U 


(1 — q> (1 2° 


rsin@ = (29) 
that solution of (28) being chosen for which rsin@ remains finite as r > a. 
Thus, the displacement of the stagnation point on the sphere is 


TPO = (30) 


2 


lim(r sin @) = 


r>a 


* Actually, r d@ becomes less small as r —- a, as the solution to be obtained shows, 
but since v, - 0 as r — a the argument is unaffected. 
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where table 1 has been used, while the displacement 4 of the stagnation 
streamline far upstream of the sphere (the main object of our investigation) 
is 


ihe I(s) ds 
6 = lim(rsin@) = | (31) 


To 
‘lo compute 6 numerically we substitute s = a cosee in (31), giving 
Ud ‘it D(a cosec ~)cosec x cot x 
Aad Aa(1 —sin®x)!* 
If Simpson’s rule with interval }7 is used in (32), we obtain 
U6 1 {s?D(s)) D(av2) v2 
— = = lim —— + 
D(a) .. COS % 
+ (a) lin. —————,, (33) 


A@ goin — 


1 Dlav2) | 2 
1 
= + 2-359 + 0-647) 
= 0:8, = 0-9. sili 


The use of a large interval }7 for Simpson’s rule may appear drastic, 
but the numerical value of the terms, not far from being in the ratio 1:4: 1, 
indicates a smoothness in the integrand which justifies the conclusion that 
Us/Aa®? = 0-9 to one place of decimals. ‘The integrations with respect 
to #, unless 19 values had been used and a few more interpolated here and 
there, were likely to have caused more error than those with respect to 
py or s which used 5 and 3 values respectively but had far greater smoothness 
in the integrands. Greater accuracy than one place of decimals is not 
required as the theory is only an asymptotic one for small Aa/U in any case. 

It is really a quadruple integration which has been performed to get 
(34). One has integrated with respect to s a function part of which consisted 
of integrations with respect to p, and @ of functions which in D were derived 
as derivatives of the ‘ drift function’ t, which is itself defined as an integral. 
The result of all this work, which it has taken three papers to expound 
(although the first two, D and I, each had, in addition, a more general 
application), is principally the numerical value (34). However, it will 
be seen in the next section that some useful qualitative results can also be 
found by studying intermediate steps of the calculation. 


4. DIscUSsION OF RESULTS AND COMPARISON WITH HALL’s THEORY 

It is helpful to begin discussing the results of §2 and §3 by comparing 
them with those of Hall’s much simpler theory of the same shear flow past 
a sphere. ‘This is based on one principal approximate assumption 
(H, p. 145): that on the plane of symmetry z = 0 the velocity gradient 
év_/ozx takes the same value as in the primary (irrotational) flow. ‘This 


j 
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assumption determines the velocity tield (v,, v,) in the plane of symmetry 
as the sum of the primary velocity field and an additional velocity field 
satisfying the two-dimensional equation of continuity 
ov ov 
+> = Q, (35) 


OX 


The only non-zero vorticity component 


— A ~ 
OX oy 


in this plane can also be determined from Helmholtz’s equation 


since the value which ov./dz takes on the plane zs = 0 has already been 
assumed. ‘lhe two-dimensional velocity field (v,,, v,,) can then be determined 
without reference to conditions outside that plane, from a knowledge of its 
divergence and curl, from the boundary condition on the body, and suitable 
boundary conditions at infinity. 

The latter boundary conditions have to be altered somewhat from their 
original form. "l'wo-dimensional flows about obstacles have the well-known 
lack of uniqueness arising from arbitrariness in the circulation; Hall avoids 
this(tacitly, in H, (2.15)) by selecting only the solutionin which the disturbance 
velocities fall off as fast as the inverse square of the distance from the origin 
(as indeed they must in the real flow); other solutions, with additional 
vorticity inside the sphere, would have disturbance velocities falling off 
as the inverse first power. On the other hand, to obtain existence of 
solutions, the condition that disturbance velocities tend to zero as r-> 
has to be applied only in a bounded interval of y (H, (3.44)), but again this is 
reasonable because the primary flow is a sensible first approximation to 
the motion only in such a region. 

With these assumptions the equations are solved by means of a power 
series in the non-dimensional shear parameter da/U. ‘The coefficient of 
Aa/U is an approximation to the secondary flow discussed in §2 and §3 
above. ‘The coefficients of (Ada/U)? and (Aa/U)%, also obtained, are 
approximations to what may be called the tertiary and quartary flows. 

In order to form an estimate of the value of Hall’s approximate 
assumptions, we may compare his results for the secondary downwash 
function D(s) with those obtained in §2 and §3. He gets, in our notation, 


D(s) 


To check this, see H, equations (2.5), (2.9), (2.11), (3.16), (3.17) and (3.21); 


in Hall’s notation, equation (38) is 


dO 
\= 


: 4 
{ U, 
(9) 
; 
| 
| 
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‘Table 2 compares values of the secondary downwash function given by 
this paper and by Hall’s theory (that is, by (38) above). It is seen that 
Hall’s method of approximation does not lead to any large departure from 
our more exact values. His downwash is an overestimate; so, therefore, 
is his inferred value of U6/Aa*, which, as deduced from (38) and (31), 
is 1:24. For comparison, if we integrate (31) approximately by means of 
equation (33), using Hall’s values of D(s) given in table 2, we obtain 
U6s/Aa*® = 1:20. Neither result exceeds our value 0-9, obtained with far 
greater labour, by as much as 40°,,. 


D(a)/Aa_ | D(ay2)/Ada} lim(s*D(s)/ Aa?) 
This paper 0-97 0-33 0-42 
H | 1-39 0-45 0-50 
Table 2. 
Now, in discussing the reasons for the success of Hall’s approach to the 


problem of estimating secondary downwash and the displacement effect, 
it must be remembered that he uses the correct secondary vorticity field in 
the plane z = 0. (For, in calculating w, to the first approximation, one 
may legitimately put ¢v./ez in (37) equal to its primary-flow value; and 
indeed Hall’s expression for w, in the plane z = 0 agrees with our equation 
(8).) His secondary flow field in that plane, however, is simply one velocity 
field which satisfies the boundary condition and has the right vorticity in 
the plane (to be precise, it is the one whose two-dimensional divergence 
vanishes). In reality, the secondary velocity field can be determined only 
from the complete secondary vorticity field, both on the plane and off, 
as in §2 above. It is easily checked that the two-dimensional divergence of 
this field in the plane z = 0 will not in general be zero. 

Thus, Hall’s assumptions are equivalent to the view that the main thing 
is to pick a velocity field which has the right vorticity locally, on the assump- 
tion that the vorticity field near the plane z = 0, with its image vorticity, 
will be much more potent in generating the velocity distribution on z = 0) 
than is the vorticity field away from z = 0 (and, in particular, the secondary 
trailing vorticity) in combination with its images. Now, this statement 
was already verified far upstream of the body in D, $3 and corrigenda, from 
which it appeared that such an approach would give a value of the secondary 
downwash only 20°, above the exact value far upstream, as is now confirmed 
by the results in table 2. For it was shown (D, (19) and (23)) that the velocity 
tield far upstream which results directly from the local vorticity distribution 
is asymptotically : 


A(V,+V,)f x 
(40) 


but that the complete asymptotic form includes also a term associated with 
the trailing vorticity, which (D, (86)) changes the V’, + V, in the downwash 


: 
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implied by (40) into V,, + 5V,, as was used in (22) above. ‘Thus in the case 
of a sphere (for which I’, = }V,) expression (40) is an overestimate by 
exactly 20°,,, which identifies it with Hall’s result in this region. 

Nearer the sphere itself Hall’s theory still gives reasonably good results 
tor the secondary downwash. ‘lo explain this it is desirable to find out 
how much of the secondary downwash calculated at s/a = 1 and v2 in §2 
comes from trailing vorticity alongside and downstream of the sphere, which 
was ignored in Hall’s theory. 

In addition, there is an independent interest in asking how much 
contribution comes from trailing vorticity in the region where a turbulent 
wake would be found in the real flow, because of course the secondary 
vorticity distribution which has been used would be seriously inaccurate 
in this region. Accordingly, the contributions to D(a) and D(av2) from 
vorticity components w,, and w,, on parts of the streamlines py/a = 0-25 
and 0-5 with @ — 90, and on parts of the streamline p,/a = 1 with 6 < 45, 
were separately evaluated. ‘Their contribution to U6/da? was found to be 
exclusively negative, but amounted only to —0-06 altogether. ‘There 
is some indication here that the applicability of the solution will not be 
much affected by the inaccuracy of the assumed distribution of vorticity in 
the wake. 

The remaining contributions to U6/ Aa? from w,,. and w,, are all positive ; 
examination shows, however, that the principal contributions come from 
fairly near the upstream axis 6 = 7. For example, positive contributions 
from the region 120 amount to + 0-06, just cancelling the (also small) 
negative contribution from the wake region. 

It is concluded that trailing vorticity cannot be regarded as responsible 
for the downwash on the axis upstream of the sphere, and this is compatible 
with the good accuracy of Hall’s approximation. Physically, the result is 
due mainly to cancellation of the effect of trailing vortices by the effect of 
their image vortices. 

Having examined the detailed computations to dispose of the suggestion 
that trailing vorticity is responsible for the displacement effect, we may 
now ask what, physically speaking, is principally responsible for it. Again, 
the computation in $2 and §3 appears to give a clear answer, namely: of 
the remaining vorticity and its images, the images make the main contribution 
to the downwash function D(s) and so to the displacement effect. 

To show this, one repeats the calculation of §2 and §3 with the terms due 
to w,, and w,, for 6 — 120° omitted (these terms were described above as 
due to trailing vorticity, and their net contribution to U6/ da? was found to be 
zero), and divides the remaining terms into those due to vorticity outside 
the sphere and those due to image vorticity. ‘Then the contribution to 
U8 Aad? of the former is found to be 0-08 and the contribution of the latter 
0-81. 

Actually, the vorticity outside the sphere makes its contribution almost 
wholly in the region far upstream; in fact the result lim{s?D(s)/Ada*} = 


12 


& 
is entirely due to vorticity outside the sphere, and this contributes 0-11 to 
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U6/Aa? in the calculation. On the other hand, downwash near the sphere 
is almost entirely due to image vorticity; the downwash at s/@ = 1 due 
to vorticity outside the sphere is only 0-15A4a (compare 0-82Aa due to image 
vorticity), and that at s/a = V2 is —0-03Aa (compare 0-36Aa due to image 
vorticity); their combined contribution to U6/Aa? is — 0-03. 

To sum up the results of this section, the expected effect of trailing 
vorticity (or more precisely of all the vorticity w,, and w,, for 6 < 120°) 
is cancelled out by that of its image vorticity; the effect of the remaining 
vorticity outside the sphere, on the other hand, is small compared with that 
of its image vorticity. 

The images of the undisturbed vorticity distribution (0, 0, — A), and of 
the ‘ring-vortex’ part w,, of the change from the undisturbed distribution, 
are especially potent in producing downwash, as is obvious geometrically 
from the rule for constructing image vortices, and they account for a part 
0-68 in U8/Aa*. images of the vortex elements w,, and w,, in @ > 120 
by contrast account for only 0-13. 


5. THe PiroT-TUBE PROBLEM DISCUSSED IN THE LIGHT OF THE RESULTS 
FOR A SPHERE 

Now, calculations like those of this paper are difficult to apply to the 
problem of the Pitot-tube displacement effect, both because of the substantial 
ditference of shape between a sphere and a Pitot tube, and because the 
calculated amount of the secondary flow, and hence also of the displacement, 
increases linearly with the shear parameter da/U; by contrast, in the range 
of Aa/U at which experiments have proved possible, the displacement 
has been found to vary little with Aa/U. 

On the first difficulty, Hall points out (H, p. 146) that, if a sphere were 
used as a Pitot tube, then, provided the round opening in the front of the 
sphere were large enough* to include the displaced position of the stagnation 
point (see (30) above for our estimate of this position), the pressure measured 
would be close to the stagnation pressure on the dividing streamline. He 
goes on to argue that this ‘spherical Pitot tube’ would in many ways be 
equivalent to an ordinary Pitot tube of somewhat smaller diameter. ‘This 
is well borne out by the subsequent work of Livesey (1956), who tested a 
Pitot tube with a hemispherical nose and a ratio d,/d, = 0-5 of the internal 
and external diameters. ‘The studies of $4 indicate that the shape of the 
front of the tube, which agrees with that of the sphere, is more important 
than that of the rear (which extends as a long cylinder instead of being 
terminated), since the trailing vorticity is unimportant, and so this tube is 
probably a lot closer than conventional ones are to the ideal ‘spherical 
Pitot tube’. (The main difference between the secondary flows in the two 
cases will be further discussed at the end of this section.) ‘The measured 
displacement was 0-10d,, as compared with 0-16d, obtained by Livesey 


* He notes that a small hole would in any case be inacceptable because of the 
resulting sensitivity to yaw. 
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(as well as by other workers) on tubes of the form used by Young & Maas 
(1936). ‘This observed difference is in agreement with Hall’s suggestion. 

Hall then goes on to meet the second difficulty by saying that the 
displacement function which we have to explain is not a constant, but rather 


is of the form 
Aa 


where sgnx is +1 when x > 0 and is —1 when x < 0. (For the dis- 
placement is in the direction y increasing if A > 0 and in the direction 
y decreasing if d <0.) In (41) the observed value of the constant C is 
about ()-2 (since the external diameter d, of the ‘spherical Pitot tube’ just 
discussed is 2a). Now, the discontinuous behaviour (41) is not plausible 
on theoretical grounds, and it is much more reasonable to suppose that 
the true behaviour is something like 


Aa 
= / 2 
( tanh( a} (42) 


say (where any odd function tending to 1 at + ~ could really be substituted 
for the tanh), but that A is large enough so that the tanh takes the value 
+1 or —1 for all da/U at which accurate measurement is possible (and 
indeed for small da/U the effect measured is so small that the experimental 
points become intolerably scattered). If (42) were correct the limit of 
U8/Aad* as Aa/U —- 0 (found in $3 to be 0-9) would be CA, so that values 
C = 0:2 and A = 4:5 would be consistent with the results of this paper. 

To test these ideas one should find the next term in the expansion of 
d/a in powers of Aa/U, which according to (42) would be 


Aa 1 Aa\3 Aa 


with the values of C and A just suggested. ‘The evaluation of the next term 
therefore requires a study of both the tertiary and quartary flows (that is, 
both the square and cube terms in the expansion of the velocity field in 
powers of Aa/U). 

Accordingly, Hall evaluates these, with the aid of his basic approximation 
(described above at the beginning of §4). After fairly lengthy calculations, 
he obtains for the displacement effect 


Aa 


‘The coettcient of (da/U)* is of the right sign, but is tive times too small 
for agreement with (43). Now, some of this disagreement may be due to 
selection of the special function tanh in (42); if one took simply a cubic 
levelled off constant beyond its maximum, the coefficient 6-1 in (43) would 
be replaced by 2:7. However, even with this form, the effect of quartary 
upwash on the displacement is still definitely underestimated by Hall’s 
theory. 


: 
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This is hardly surprising. In the first place, the cumulative error in 
applying Hall’s approximation three times to obtain successively the 
secondary, tertiary and quartary flows must greatly exceed the error in 
applying it once only, which already produced an error in the displacement 
of the order of 30°. Again, Hall’s approximation gives erroneously only 
the divergence, not the curl, of the secondary flow field in the plane of 
symmetry; but both divergence and curl will be erroneous in the case of the 
tertiary and quartary flows. 

A more clear-cut argument, which also appears to explain the sign of 
the error in Hall’s calculation, is provided by some discussion in D, §3, 
which is greatly amplified and extended in a forthcoming paper (Lighthill 
1957 b). 

There it is shown that far upstream the tertiary flow tends to zero more 
slowly than the secondary flow, and that similarly the quartary flow is greater 
in magnitude than the tertiary (not, in fact, tending to zero at all). Hall’s 
equations give a smaller order of magnitude for the quartary upwash, at 
least far upstream, than that indicated by the exact theory. This indicates 
that the true reduction of the displacement by the quartary flow may be 
greater than Hall predicts. 

Actually, the successive-approximation sequence is not uniformly valid 
for large ry. In this region the disturbances can best be treated as a small 
perturbation of the exact parallel flow, which leads to their expression as 
a Hankel transform of suitable solutions of the steady inviscid case of the 
Orr—Sommerfeld equation for the shear layer. 

The solution so obtained can be represented as the sum of secondary, 
tertiary, quartary flows and so on, with downwashes of orders s~*, s~}, 1, 
respectively, only for moderately small s. For larger s the downwash 
behaves as a more complicated function, which ultimately falls off like 
sass «. The solution of this paper which assumes that the falling-off 
is like s-? all the way therefore overestimates the displacement. It can be 
shown (Lighthill 1957 b) that an approximate form of the correction to 
5 due to departure from the secondary flow for large s is 


(45) 
where the shear A has been supposed a function of y (measured from an 
origin at the centre of the sphere). ‘lhis term of order a? is negative, and is 
intermediate in order between the secondary flow term 0-9Aa?/U of this 
paper and Hall’s quartary-flow term of order a*. ‘This again gives a reason 
for the fairly rapid turning over of the graph of 5/a against Aa?/U. 

There is one more difficulty in the comparison between the theory for 
the sphere and the experimental results, which can be investigated effectively 
in the light of the theory just described. It can take two forms. 

First, if an actual sphere were used as a Pitot tube, the fact that the sphere 
experiences drag, and has the wake associated with its drag, modifies the 
upstream character of the primary irrotational flow, especially far from the 
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body. ‘There it is asymptotic not to a ‘doublet’ motion but to that due 
toa source of strength D/pU, where Disthe drag. Accordingly the upstream 
behaviour assumed in our solution is incorrect. ‘The weight of this objection 
should not, however, be overestimated. For a sphere of radius a the doublet 
strength is 27la*, while the source strength even under the high-drag 
conditions obtaining with laminar separation (with, say, Cp = }) is only 
\7Ua*. Hence the source produces a velocity reduction upstream of the 
body greater than that due to the doublet for r > 16a only, so that it could 
hardly be expected to be important. 

Secondly, for a real Pitot-tube form, the upstream behaviour of the 
primary irrotational flow will again be asymptotically that of a source, this 
time of strength 7Ua®, where a is the external radius. Obviously it is more 
questionable in this case whether the results calculated for the sphere are 
applicable, even though the Pitot tube may have a hemispherical nose, 
because of the larger source strength. In both cases the problem is made 
more serious by the result (D, (29)), which gives 


Am 


D(s) (46) 
where m is the source strength. If the downwash really satisfied (46), 
then the integral (31) for the displacement would be logarithmically infinite. 

However, it is shown in the forthcoming paper already mentioned that, 
when the solution for large s is replaced by a uniformly valid approximation, 
a finite value for 6 emerges. It is shown that the displacement so obtained 
is the same as that deduced from the secondary flow alone if, in the latter, 
the integration be carried only to the value s = s,, described as the secondary- 
How cut-off. In other words the displacement is the same as if the secondary 
flow penetrated only a distance s, ahead of the source. ‘The distance s, 
is of the order of the width of the shear layer; an approximation to it is 
given by 


(47) 
Og s,. = | — —sgny dy, 7 
which makes s, in a sense a‘ geometric mean’ of the difference in y-coordinate 
between points in the shear layer and the source. 

‘These considerations lead to a term 


Am ope dA y 


in the expansion of 8 in powers of a, where in (48) the value of m appropriate 
to a solid Pitot tube of external diameter d, = 2a has been inserted. ‘The 
term (48) of order aloga in 6 a is in addition to the previously found term 
of order a. ‘The logarithmic term certainly is a non-negligible addition in 
this case; when the integral is as much as 24, it would increase the coefficient 
of da/Ubyabouta quarter. Inthe other case, when only the source resulting 
from the wake drag is present, the effect is not very important. In both cases, 
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however, the effect is not a complete change of order of magnitude but 
rather a mere quantitative increase. 

‘To sum up, we may state as the main conclusion of this section that the 
original Young & Maas suggestion of a discontinuous functional dependence 
(41) of 6/a on Aa/U is without theoretical support. All the theories agree 
in predicting a continuous functional dependence, which one can only 
assume (in the light of experimental results) becomes slowly-varying for 
values of Aa/U in excess of a fairly modest limit, an assumption which in 
the light of the full discussion is seen to be not inconsistent with the theories. 

We may ask, in conclusion, whether the work of this paper gives any 
explanation of two results which were noted as especially interesting in the 
introduction: the great reduction of the displacement effect with supersonic 
flow (Johannesen & Mair 1952) or with sharp-lipped tubes (Livesey 1956). 
Certainly no precise answer is possible in either case, but it is probably 
relevant to both facts that, as shown in $4, the image vorticity is the main 
contributor to the diplacement effect. 

Thus, in supersonic flow about a Pitot tube, the extent of regions of 
vorticity whose image vorticity could affect the downwash on the dividing 
streamline is much reduced—partly because the region of subsonic flow 
is cut off by a shock wave a short distance upstream of the orifice, and partly 
because the flow on streamlines as far away from the axis as py/a = 1 (which 
were found in the calculations of §2 to make substantial contributions to 
the downwash) becomes supersonic very soon, and ceases to have a domain 
of influence which includes the dividing streamline. 

Similarly, one may imagine that the ideal Pitot tube with internal and 
external diameters equal (to which Livesey’s sharp-lipped tube is a good 
approximation) produces very little displacement because image vorticity 
in sucha tube is reduced, both in magnitude and in effectiveness for producing 
displacement. For, first, the undisturbed vorticity distribution (0, 0, — A) 
has no image vorticity at all in this case. In other words, the ‘first part’ 
of the secondary flow, as defined in D, §3 and in §2 above (and which for 
a sphere has the potential (3)), vanishes, since the uniform shearing motion 
by itself possesses no velocity component normal to the surface. 

As to the image system of the vorticity change @,, one may note that this 
will include a uniform vorticity (0, 0, +A) from some point of the tube 
downstream inside the tube; for in this region all motion, and hence also 
all vorticity, is doubtless absent. However, the effect of this vorticity 
trapped inside the tube is negligible outside it; even in the corresponding 
two-dimensional case it falls off exponentially, by a factor of e* = 23 in one 
diameter. 

Thus, we are led finally to ask about the effect of such disturbance vorticity 
as is placed similarly to that which produced the displacement in the case 
ofasphere. Asin that case, the effect of the trailing vorticity will be cancelled 
out by the effect of its images; and the vorticity upstream of the Pitot tube 
may be expected to make its effect largely through its image vorticity. 
However, since there is no room inside the surface in the case of a shape like 
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this, the image vorticity must be rather remote—in mathematical language, 
it is on another sheet of a Riemann surface (of three dimensions). 
Accordingly, its image effect upstream may be expected to be less than in 
the case of a more full-bodied shape. 

Again, we have not explained the results by the theory, but merely shown 
that the theory is not inconsistent with them. Clearly, much more work 
will be required before theories of shear flow become fully satisfactory and 
widely applicable. However, the attempt to achieve this end seems to be 
worth making. 
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On the existence of higher than normal detonation 
pressures 


By DONALD R. WHITE 


General Electric Company Research Laboratory, Schenectady, New York 


(Received 25 March 1957) 


With the increasing popularity of combustion driving in strong shock 
tube experiments, the possibility of detonation must be considered seriously. 
In addition to the extreme and relatively well-known conditions associated 
with a detonation, there exists a phase during the transition from a flame 
to a detonation during which pressures and velocities of propagation can 
exist which are higher than those associated with the fully formed detonation. 
This has been discussed theoretically by Oppenheim (1952) and has been 
shown experimentally by several workers. A luminous front having a higher 
than detonation velocity has been observed by Bone, Fraser & Wheeler (1936) 
and others. J. B. Smith (1949), using calibrated burst diaphragms at the 
end of a pipe, has observed reflected pressures in fuel gas-air, acetylene—air, 
and hydrogen-air mixtures which were approximately four times greater 
than those associated with a reflected normal detonation. ‘Turin & Huebler 
(1951), using quartz pressure transducers in the side of a tube, observed 
pressures during this transition process about three and a half times higher 
than those at a later time when detonation was fully developed. ‘They 
worked with ‘Toledo natural gas. Mooradian & Gordon (1951) also have 
clear indications of this phenomenon with hydrogen—air mixtures. 


‘Tube ‘Time to half 

Peak pressure ; Impulse 
length maximum pressure 
24 ft. no detonation a _ 
52 ft. 2600 Ib. 200 microsec 0-8 Ib. sec/in.* 
47 ft. 710 Ib. in.* 1000 microsec 1-0 Ib. sec in.” 


Table 1. 


Several experiments have been run in this laboratory using a quartz 
pressure transducer (SLM model PZ 14) mounted in the end flange of a 
3} in. square tube. A stoichiometric hydrogen-air mixture at atmospheric 
pressure was ignited near one end of the tube with a spark plug, and the 
pressure on the other end of the tube was measured as a function of time. 
This pressure record was then integrated between the time of impact and 
the time the pressure dropped to 150 Ib./in.* abs. to obtain the impulse 
applied to the end of the tube during this interval. ‘lhe pressure level 
corresponding to constant volume combustion is about 120 Ib./in.? abs., 
and that due to a reflected normal detonation wave is on the order of 500 
or 600 Ib./in.?.. Our measurements are given in table 1. 
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The 32 ft. length was marginal with respect to occurrence of this over- 
detonation when the ignition was at the extreme end. However, moving 
the spark plug 32 in. from the end resulted in this over-detonation on each 
run at this length. Over-detonation was accompanied by a much more 
substantial noise than was detonation, like the ditference between blows 
by a sledge hammer and a carpenter’s hammer. 

This over-detonation phenomenon may have escaped the attention of 
some experimenters using combustion driving in shock tube experiments. 
It must be noted that not only are the limits for the occurrence of detonation 
rather uncertain, but also that the first departure from the expected constant 
volume combustion may result in a pressure loading much more severe 
than that which would result from a fully developed detonation. ‘The long 
distances noted above do not offer comfort since the use of oxygen instead 
of air greatly reduces the formation distance for detonation. ‘This hazard 
may be reduced by using a number of points of ignition and by choosing 
a gas loading sequence such that the mixture does not go through a condition 
which is more readily detonable than the final mixture desired (i.e. loading 


of hydrogen last). 
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Viscous Flow Theory, I—Laminar Flow, by Suin-I Par. Princeton: 
D. Van Nostrand Co., 1956. 384 pp. $7.75 or 58s. 


In the study of real fluids, this has been the half-century of boundary 
layer theory. It is remarkable how Prandtl’s discovery still dominates 
research in viscous flow, as it has done since 1904, though one may suspect 
that the possible extensions of boundary layer theory are at last nearly 
exhausted, and that the really exciting discoveries await us in problems in 
which the boundary layer, if present at all, is but one component of the 
flow field. Oseen swam against this current for years, not finding it 
necessary in his Hydrodynamtk of 1927 even to mention Prandtl or his 
Grenzschicht. But with a few other exceptions the history of viscous 
flow theory since the turn of the century has been in large part the record 
of the development of boundary layer theory, following paths that to a 
remarkable extent were at least sketched in outline by Prandtl and his 
students in the first few pioneering years. 

It is perhaps not surprising, then, that the book under review, which 
is the first of a two-part work on viscous flow theory (the second being 
devoted to turbulent flow), is concerned mainly with the theory of the 
laminar boundary layer. Indeed, aside from two short and rather dis- 
appointing chapters on other topics, it is devoted entirely to a thoroughly 
modern survey of laminar boundary layer theory, together with the 
requisite introductory chapters. It is not intended to be a comprehensive 
treatise in the style of Schlichting’s Grenzschicht-Theorie, but a textbook 
for advanced students of engineering and aerodynamics. ‘The account is 
purely theoretical, with no enlightenment ever sought from experiment. 

The great virtue of the book is that it is the first to assign a central 
role to compressibility and heat conduction. In this respect it is strongest 
where Schlichting’s treatise is weakest; for, even in the expanded English 
translation, Schlichting considers an incompressible fluid from the outset, 
and defers to the advent of high-speed flight only in two appended chapters. 
(Goldstein’s Modern Developments in Fluid Dynamics, written as it was at 
the dawn of the age of compressibility, could not reasonably have attempted 
a unified treatment of compressibility and heat conduction.) ‘lo be sure, 
Dr. Pai is forced more often than not to retrench in view of mathematical 
complexity, and to consider a particular problem only in the special case of 
incompressible flow. But his eye is constantly on the compressible fluid. 
This is undoubtedly the modern approach, to be taken by future writers. 

If one thus breaks with tradition and starts from the Navier-Stokes 
equations for a compressible heat-conducting gas, it is imperative that he 
understand precisely how compressibility enters the problem, and how 
the often necessary specialization to incompressible flow is to be made. 
The supersonic aerodynamicist, who automatically associates compressi- 
bility with Mach number, may be surprised to realize that for a heat- 
conducting fluid there is a second way in which compressibility can appear. 
No matter how slight the Mach number, heating of a body will induce a 
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temperature ficld in the flow and consequent variations of density. 
Unfortunately, such important distinctions are not always clearly made in 
the present book. Nor are they perhaps entirely clear to the author, who 
remarks (p. +9): ‘* Another current practice in finding the temperature 
distribution approximately is to drop the viscous dissipation terms in the 
energy equation, and then to solve for the temperature distribution. 
Strictly speaking, such an approximation is not logical, because the terms 
of heat conduction and vicous dissipation are of the same order ot 
magnitude. But for practical purposes, when both the velocity of the 
flow and the temperature gradients are small, such an approximation 
seems to give reasonably good results.” 

Now the fact of the matter is that conduction and dissipation are not 
of the same order of magnitude in general, and that there are important 
classes of problems in which it is an entirely logical approximation to drop 
the dissipation term (and the more so the larger the temperature gradients). 
Consider for definiteness a stream of perfect gas flowing at Mach number 
M past a heated body. Introducing dimensionless variables in the energy 
equation (or solving an example) shows that both the dissipation and the 
enthalpy changes due to pressure gradients are proportional to M?*. 
Hence dissipation and also pressure terms in the energy equation can be 
neglected strictly at J/ = 0 or approximately when JJ? is very small. If, 
in addition, the maximum temperature difference AT is small compared 
with the temperature 7 itself, then density variations are slight, and to a 
first approximation the momentum and energy equations are not coupled, 
so that the heat is transported by the incompressible velocity field. (‘This 
is readily verified by expanding all flow quantities in powers of AT T and 
retaining only linear terms.) 

If the body is insulated (the thermometer problem), viscous dissipation 
is the only source of heat, so that a problem exists only for W> 0. If VP 
is small, the heat 1s again carried by the incompressible velocity field to a 
first approximation. ‘The combined problem in which heat is supplied 
both by transfer from the body and by dissipation can be handled with a 
double series expansion, which shows that to a first approximation for 
small \7 7 and .W> the solution of the thermometer problem can simply 
be added to the solution of the heat transfer problem for WV = 0, and that 
non-linearities as well as coupling between the momentum and energy 
equations arise only in the second-order terms. 

It should be pointed out that the above conclusion that dissipation is 
strictly negligible at 17 = 0 holds only for flows with a definite stream. In 
problems of free convection the dissipation can be neglected only as a 
first approximation for slight heating, because the velocities are pro- 
portional to A7’ 7 and the dissipation is consequently of order (AT T)?. 
It must be admitted that these matters are not well explained in most 
existing books on viscous flow theory, or in many research papers. (‘They 
are very carefully discussed in Lagerstrom’s forthcoming article on laminar 
flow in the Princeton series, and briefly but precisely in Lighthill’s paper in 
Proc. Roy. Soc. A, vol. 202, 1950, p. 359.) 
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The first hundred pages of the present book are introductory in nature, 
covering basic thermodynamics and kinetic theory, the laws of friction and 
heat conduction, the derivation of the equations of motion, a number of 
illustrative exact solutions (mainly incompressible), dynamic similitude, 
and a general discussion of the Navier-Stokes equations, including Burgers’ 
model, linearization, and the limits of low and high Reynolds number. 
‘There follows a short chapter on slow motion, which has a curiously old- 
fashioned tone. Indeed, there are no references to work done since 1910, 
which is a pity, since the realm of Stokes and Oseen flow is one in which 
some interesting discoveries have recently been made, with more expected 
to follow. 

The main bulk of the book, some 200 pages, is devoted to a generally 
admirable survey of laminar boundary layer theory. ‘This is the field of 
the author’s own research, is clearly where his heart lies, and is by far the 
strongest part of the book. ‘lhe emphasis is on compressible flow where 
possible, including such recent matters as the hypersonic boundary layer 
and interaction with a shock wave. Successive chapters treat the boundary 
layer equations, exact solutions, approximate methods, axisymmetric and 
three-dimensional boundary layers, unsteady boundary layers (including 
a sketch of the stability theory), and boundary layers with suction and 
injection. Aside from Kuerti’s survey in the second volume of Advances 
in Applied Mechanics, there is no comparable summary of compressible 
laminar boundary layer theory, and I would recommend it with only minor 
quibbles. I would only protest that more than the necessary two boundary 
conditions are provided for the energy equation on pages 148-149, and that 
of these the requirement of equal normal temperature gradients in the gas 
and body at their juncture is incorrect. (Note also the unfortunate typo- 
graphic loss, on page 178, of a radical sign in Blasius’ famous formula for 
flat-plate skin friction.) In particular, the author displays considerable 
skill in summarizing the various complicated extensions to a compressible 
fluid of the Karman-Pohlhausen approximation, 

He breaks away from the boundary layer approximation again in the 
last chapter, which is entitled ‘* Linearized theory of viscous compressible 
fluid”. Although mainly a survey of Lagerstrom, Cole, and ‘Trilling’s 
study of compressible Oseen flow, it also includes, among other incidentals, 
an overly-condensed account of Kuo’s second approximation to Blasius’ 
solution for the boundary layer on a finite flat plate, from which the main 
result—the numerical coefficient of the second term (it is a, = 4:12)—is 
unaccountably omitted. Unfortunately, the author has thought it necessary 
to copy some ten pages of his book literally word for word from the report 
of Lagerstrom et al. (In his only major departure from that text, he 
thoroughly garbles the characteristic condition for a second-order partial 
differential equation.) Now it is perhaps a matter of taste whether one 
chooses to expound another’s theory by quoting him verbatim and at length; 
this seems to be a growing trend in the literature of fluid mechanics, for I 
have found three well-known writers doing so within the last year. But it 
is certainly another, and quite indefensible, matter when the present author 
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concocts introductory pages that are ostensibly his own by assembling a 
pastiche of sentences and paragraphs lifted from the Lagerstrom report. 
‘This is, then, a book with some serious flaws, written in a style that is 
often derivative, and concerned more with examples than ideas. Despite 
its title, it concentrates on the boundary layer approximation to the virtual 
exclusion of other facets of viscous flow theory. Yet, admitting these 
limitations, it is the best available survey of laminar boundary layer theory 
from the modern point of view. ‘The student who arms himself also with 
at least the first of Goldstein’s two masterly volumes, to provide the 
physical insight and appeal to experiment that the present book lacks, 
should find it a useful guide to that subject. M. D. Van DyKE 


Molecular Flow of Gases, by G. N. Patrerson. New York: John 
Wiley & Sons, 1956. 217 pp. $7.50 or 60s. 


The molecular viewpoint in studies of gas flows is by no means a new 
area of interest. It goes back at least to the studies of Kundt, Warburg 
and Maxwell (1875-1879) on low-speed slip flow, and the later work of 
Knudsen and Smoluchowski (1910-1911) on low-speed free molecule flow. 
However, the increase in the amount of research in molecular gas flows 
since World War [I] is most impressive. Over and above a proportionate 
share of the post-war growth of all branches of science, the increasing 
importance of the molecular viewpoint in many problems has attracted the 
attention of more and more scientists from the disciplines of physics, 
chemistry and astronomy. 

There is a wide range of problems requiring a consideration of the 
molecular structure of moving gases. Such considerations are necessary 
in the study of colliding galaxies, supernovae, stellar atmospheres, and 
meteors, in the development of high-speed, long-range missiles, in the 
development of controlled thermonuclear reactions, and in the study of 
shock-induced chemical reactions. An illustration of the kind of problem 
to be dealt with is provided by a simple shock transition. For weak shock 
waves in a moderately dense monatomic gas, the shock transition zone 
(the region between the shock front and the attainment of equilibrium in 
the shocked state) is so thin as to defy detailed experimental study, the gas 
behaves ideally, and the transition can be dealt with adequately by regarding 
it as a discontinuity in a continuum fluid flow problem. Conservation 
laws are here sufficient for a determination of the change in flow properties 
across the transition. At low densities in a monatomic gas the density 
profile through the shock thickens in proportion to the increased mean 
free path in the gas, for a certain number of collisions are required for 
equilibration of translational energy of the molecules. (It is interesting 
to note here that a shock profile derived from the Navier-Stokes equation 
with the continuum viewpoint represents the observed transitions better 
than the best solutions of the Boltzmann equation yet available.) With 
stronger shocks in a diatomic gas the transition zone is thickened further, 
as the gas requires many more collisions to equilibrate rotational and 
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vibrational degrees of freedom than to equilibrate translational kinetic 
energy. Here we must not only consider the Boltzmann equation with 
inelastic collisions, but must modify the ideal gas form of the internal 
energy to include thermodynamic properties of the specific gas under 
study, in order to specify correctly the final equilibrium state. For still 
stronger shocks we must include molecular dissociation, ionization and 
other gas phase reactions. Much work needs to be done before we 
understand on a theoretical basis even such a basic problem as the shock 
transition in molecular gas flows. 

There appear to be several stages of sophistication involved in 
considering gas flows on a molecular basis. ‘There is firstly a purely kinetic 
viewpoint in which one assumes elastic collisions and a model for the 
intermolecular forces in the gas. If the collisions cause at most a small 
deviation from a Maxwell velocity distribution, then one can use the 
Boltzmann equation and the Maxwell equations of transfer to develop the 
fluid equations of motion, including transport properties. Boundary 
conditions at a surface in the flow also may be developed on this basis. 
A more sophisticated development would describe gas flows under the 
condition that the molecular collisions are inelastic. Here the collisions 
not only alter materially the velocity distribution function, but involve 
energy transfer into internal excitation, dissociation, ionization, recom- 
bination, and so on. Ina still more sophisticated analysis one might attempt 
(although whether it would be profitable to do so is debatable) to deal with 
turbulent flows, in which one might combine consideration of microscopic, 
random molecular motions with the ordered motion in eddies having 
dimensions of many gas mean free paths. Other more complex flow 
problems include the motion of highly ionized gases involving collective 
behaviour due to long range forces, and the motion of very dense gases in 
which one may not be able to limit the analysis to bimolecular collisions 
only and in which collisions may not be isentropic because of the finer 
details of intermolecular forces. 

Professor Patterson’s book is concerned primarily with the first stage 
of sophistication mentioned above. ‘lhe book has as its stated purpose 
the provision of a transition from the continuum to the molecular viewpoint 
in fluid mechanics, and provides perhaps as elementary an introduction to 
the kinetic theory of non-uniform gases as can reasonably be expected for 
those interested principally in fluid mechanics. It might be read with 
profit by those who would later undertake study of such works as Mathematical 
Theory of Non- Uniform Gases by Chapman and Cowling or Molecular Theory 
of Gases and Liquids by Hirschfelder, Curtiss and Bird. 

The book begins with a brief discussion of basic concepts in kinetic 
theory, including an introduction to the Boltzmann integro-ditferential 
equation and the Maxwell transfer equation. A very clear discussion is 
given of molecular collisions. Next the author describes the solution 
of these equations for gases represented by perfectly elastic spherically 
symmetric particles, having no internal degrees of freedom. Collisions 
of such particles do not alter the velocity distribution function, the 
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gas is in a steady state, and the Boltzmann equation leads directly to 
the Maxwellian velocity distribution. ‘lhe transfer equation then leads 
to the equations of isentropic tlow. Properties of a Maxwellian gas are 
then discussed. As part of the isentropic flow discussion the author 
develops at some length the non-stationary expansion wave. 

‘The major part of the book is devoted to slightly non-isentropic flows, 
i.c. flows in which small gradients in composition, temperature, and density 
occur, and in which the velocity distribution function deviates slightly 
from the Maxwellian. ‘lo describe such flows one assumes a repulsive 
force fieid for a point particle having no internal degrees of freedom. 
Using a method due to Grad, the author obtains a perturbation solution 
to first order of the Boltzmann equation for a gas not in a shady state. 
‘The resulting transport properties are similar to the Chapman-Enskog 
result. ‘The theory is then applied to weak shock transitions and boundary 
layer flows. Comparisons are made with experimental studies of slightly 
non-isentropic flows. A very brief discussion is given of the flow alterations 
resulting from real gas ettects, i.e. internal excitation of atoms and molecules, 
dissociation, 1onization, etc. Finally there is a description of the present 
state of the theory of slip flow and of free molecule flows. Appendices 
include some mathematical relations which occur in kinetic theory problems, 
a discussion of first-order ordinary and partial differential equations, and 
a concise statement of the procedure leading to the Chapman-Enskog 
result and the Burnett approximation. 

Several features and flaws of the book invite comment. ‘The author 
has achieved his goal of giving an introduction to the flow of gases with 
molecular structure. ‘The discussion of molecular collisions and_ the 
development of the flow of raritied gases are presented with particular 
clarity. ‘The exposition of expansion waves and characteristics is un- 
necessarily long, not too pertinent to the main thesis of the book, and the 
material covered is readily available to the reader elsewhere. ‘The discussion 
of gas flows in which molecules are excited, dissociated or ionized, or in 
which chemical reactions may occur, is much too brief. What discussion 
is given Is cursory and fails to convey to the reader what can be accomplished 
by an analysis of high temperature and high velocity flows in terms of 
thermodynamics and chemical kinetics. For example, one finds the familiar 
equilibrium constant applied to a dissociation reaction referred to as 
‘a characteristic quantity relevant to dissociation”’. ‘The appendices could 
be eliminated without any particular loss to the book, except that the brief 
mention of the Burnett approximation might profitably be expanded and 
incorporated into the text as part of a discussion of the general problem of 
dealing with the Boltzmann equation. ‘The author does not even hint 
that there is at present considerable interest in the molecular description 
of the flow of conducting gases in magnetic fields. 

The flaws mentioned nevertheless are minor, and do not detract 
seriously from the general worth of the book as a useful text for introduction 
of the student of gas dynamics to the molecular viewpoint. 


R. A. ALPHER 
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